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PREFACE 


This Second Annual Status Report covers the effort of the base program 
on the NASA HOST program titled "Constitutive Modeling for Isotropic Mate- 
rials" conducted under Contract NAS3-23925. The NASA program manager for this 
project is Mr. Albert Kaufman. The program manager at Southwest Research 
Institute is Dr. Ulric Lindholm. Contributors to this report are Dr. Kwai 
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1.0 INTRODUCTION 

The objective of the present program is to develop a unified constitu- 
tive model for finite-element structural analysis of turbine engine hot sec- 
tion components. This effort constitutes a different approach for non-linear 
finite-element computer codes which have heretofore been based on classical 
Inelastic methods. The unified constitutive theory to be developed will avoid 
the simplifying assumptions of classical theory and should more accurately re- 
present the behavior of superalloy materials under cyclic loading conditions 
and high temperature environments. This class of constitutive theory is char- 
acterized by the use of kinetic equations and Internal variables with appro- 
priate evolution equations for treating all aspects of inelastic deformation 
including plasticity, creep and stress relaxaMon. Model development is 
directed toward Isotropic, cast nickel-base alloys used for air-cooled turbine 
blades and vanes. Recent studies have shown that this approach is oarticular- 
ly suited for determining the cyclic behavior of superalloy type blade and 
vane materials and is entirely compatible with three dimensional Inelastic 
finite-element formulations. More efficient and accurate inelastic analysis 
of hot section components— turbine blades, turbine vanes, combustor liners and 
seals— fabricated from "age hardenable" Isotropic superalloy materials will be 
realized as the result of these developments. 

The program Is being conducted In two phases. A basic program (Tasks 
A through I) and an optional follow-on program (Tasks J through M). In the 
basic program of twenty-six months' duration, a unified constitutive model 
will be developed for the prediction of the structural response of Isotropic 
materials for temperatures and strain range characteristics of cooled turbine 


vanes in advanced gas turbine engines (Task Al. A data base of uniaxial and 
multiaxial material properties required for the constitutive model development 
will be obtained for the base material (Tasks C and E). The constitutive 
model will then be incorporated into a finite-element computer code (Task D). 
An evaluation will be made of the capability of the analytical method to pre- 
dict the structural response for multiaxial stress states (Task E) anJ noniso- 
thermal conditions by conducting thermomechanical loading and benchmark notch 
verification experiments and analysis (Task F). As a final evaluation of t^e 
analytical model, a structural analysis will be performed for a hot section 
component fabricated from the base material for simulated engine operating 
conditions (Task G). In the optional program, material property test proce- 
dures will be further developed to minimize the amount of testing required, 
and to study the possibility for estimating the material model constants from 
conventional property data (Task J). Further development of the model will be 
undertaken to consider thermal history effects and to correct any deficiencies 
indicated in the model or the computational algorithms in the code (Task K). 
In addition, the constitutive model development will be verified for an alter- 
nate material (Task L). 

The work under this program is being conducted as a joint effort be- 
tween Pratt & Whitney Aircraft (PWA) and Southwest Research Institute (SwRI) 
with technical assistance from Prof. Sol R. Bodner and Or. Kevin Walker in the 
area of constitutive model development. The work and data base generated is 
being coordinated closely with another NASA funded HOST program at PWA (NAS3- 
23288) to develop advanced life prediction techniques for isotropic superalloy 
vane and blade materials. 



I 


This report documents the efforts in the base program (Task A through 


I). Although efforts in certain tasks were completed and reported last year 
(1), some of those results will also be included in this report for the sake 
of clarity and completeness. 
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2.0 MATERIAL AND CONSTITUTIVE MODELS SELECTION 


2.1 Material Selection 

The model material selected for study in the base program is B1900+Hf 
and MAR-M247 has been chosen for the optional tasks. Both materials are ex- 
tensively utilized in the gas turbine industry and considerable benefit will 
derive from their characterization. A major factor in the selection of 
B1900+Hf is the availability of cyclic and monotonic constitutive data from a 
concurrent NASA HOST program, NAS3-23288, conducted by PWA to develop life 
prediction analysis methodology for isotropic blades and vanes. Selection of 
the same material, material processing, and specimen configurations in both 
efforts significantly enhances the data base available for life prediction and 
constitutive model development. 

The B1900+Hf material was part of a single heat, designated W-0098, 
obtained from Certified Alloy Products Inc., Long Beach, California. The 
chemical composition of this heat is compared to nominal specification in 
Table 2.1. Casting configurations, pour and mold temperatures for the 
B1900+Hf specimens were selected to assure that grain size, material structure 
and integrity match, as closely as possible, those in the PWA test specimens. 
Goals for this phase of the fabrication effort were: 

1. A grain size of ASTM No. 1 to 2 in the gage section 
when measured using standard procedures. 

2. A porosity-free casting. 

A y ' size of .9 urn in the fully heat treated condi- 


3. 


TABLE 2.1 


CHEMICAL COMPOSITION OF B190O+Kf (HEAT W-0098) 


Element 

Nominal (%) 

Heat W-0098 

C 

0.11 

0.09 

Cr 

8.0 

7.72 

Co 

10.0 

9.91 

Mo 

6.0 

5.97 

A1 

6.0 

6.07 

Ta 

4.3 

4.21 

Ti 

1.0 

0.99 

B 

0.015 

0.016 

Zr 

0.08 

0.04 

Fe 

0.35* 

0.17 

M 

0.1* 

0.04 

Cb 

0.1* 

0.08 

Bi 

0.5 ppm 

0.1 

Pb 

10.0 ppm 

0.1 

Hf 

1.5 

1.19 

Ni 

Remainder 

Remainder 


* Maximum 
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Following casting, all B1900+Hf test bars were fully heat treated ac- 
cording to the following schedule: 

Solution: 1079+14 o C (1975±25°F) for 4 hours; air cool 
Precipitation: 899+14*0 (1695±25°F) for 10 hours; air cool. 

Figure 2.1 shows the grain size of the as-cast and the fully heat- 
treated material. The y size is approximately .6 pm 7 n the as-cast material 
and is .9 pm in the fully heat-treated condition, as shown in Figure 2.2. 

2.2 Constitutive Model Selection 

A literature survey has been conducted to assess the state-of-the-art 
of time-temperature dependent elastic-viscoplastic constitutive theories which 
are based on the unified approach. This class of constitutive theories is 
characterized by the use of kinetic equations and internal variables with ap- 
propriate evolution equations for treating all aspects of inelastic deforma- 
tion including plasticity, creep, and stress relaxation. As reported earlier 
[1), the review identified more than ten such unified theories which are shown 
to satisfy the uniqueness and stability criteria imposed by Drucker's postu- 
late and Ponter's inequalities. These theories are compared on the basis of 
the types of flow law, kinetic equation, evolution equation of the internal 
variables, and treatment of temperature dependence. 

As a result of the literature survey and based upon SwRI recommenda- 
tions, the models of Bodner and Partom |2] and of Walker [3] have been selec- 
ted and approved for further study. These two models are both representative 
of the class of unified models considered in the review process but differ 
significantly in the choice of particular functional forms for the basic flow 
law, the kinetic relationship, the parameter used as a measure of hardening 
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Fully Heat Treated y* = 0.9 /in 10 P m 


FIGURE 2.2. GAMMA PRIME (y‘) SHOWS INCREASE WITH HEAT TREATMENT. 
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and the evolution equations for the Internal variables describing work harden- 
ing. Thus, a direct comparison between these two models and the experiments 
should illustrate well the consequences of a wide range in constitutive 
modeling approach. It Is also significant that both models have already found 
significant application to analysis of gas turbine materials and to hot sec- 
tion components. Therefore, they are further along in their development and 
evaluation than most of the other comparable models. 

As a result of the literature survey and of direct experience in work- 
ing with existing models, modifications to the Bodner-Partom and Walker models 
have been made to account for additional hardening resulting from multiaxial, 
non-proportional load or strain histories. Modifications have also been made 
to the kinetic equation and the evolution equations of the internal variables 
of the Walker model. The updated Bodner-Partom and Walker models will be dis- 
cussed individually in Section 3.0. 
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3.0 BRIEF REVIEW OF CONSTITUTIVE MODELS 
3.1 Bodner and Partom Model 

The elastic-viscoplastic theory of Bodner and Partom, published in 
1968 [4], was probably the first "unified" set of constitutive equations with- 
out a yield criterion or loading/unloading conditions to be developed. Those 
equations include certain physical concepts provided by the work on disloca- 
tion dynamics during the 1950 ' s and early 1960 ' s . The equations are placed in 
the context of multi-dimensional continuum mechanics which makes them capable 
of solving problems by analytical and numerical methods. One of the initial 
papers, in 1972, considered large deformations [5] and another the same year 
included isotropic hardening (6). At the present stage of development, the 
constitutive theory includes isotropic and directional hardening, thermal re- 
covery of hardening, general temperature dependence of plastic flow, and iso- 
tropic and directional damage development [7-9]. 

A summary of the current Bodner-Partom model without the damage term 
is shown in Table 3.1. As listed in Table 3.2, there are 13 material con- 
stants in the model, most of which are temperature independent, and generally 
fewer (about 9) are needed in an actual material representation. The basic 
equation for plastic strain rate as a function of stress and history dependent 
Internal variables is taken to be of exponential form. The initial work [4], 
before 1970, was based on a power law relation which seems to have been 
adopted by other investigators. Properties of the exponential function used 
in the B-P formulation that may be especially useful are its very low value, 
almost zero, at low stress levels and Its limiting value of plastic strain 
rate. As a consequence, that particular function seems to be suitable to 
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TABLE 3.1 

A SUMMARY OF BODNER-PARTOM'S MODEL 


BODNER-PARTOM'S MODEL 


Flow Law: 


- i* * i.? 


6 1j 


u 


*1j * x S ij ; *kk * 0 

with $ij • * y <*i j®Kk 


Kinetic Equation: 

°Z * °o 2 exp 
with Z » Z 1 + Z° 

°2 P * 1 2 1j *1J 





J 2 * I S ij S 1j 
X 2 - 0 2 p /J 2 

Evolution Equations of Internal Variables: 
a. Isotropic Hardening 

.1 


2 1 - Z 9 

-ir* 


where 


Z . mi [z i+ aZ 3 -Z I jw p -A 1 Z 1 

i a 1 * 2 (ot] - a) lip Sin0 
9 * cos* 1 (V^V.j) or 9 » cos' 1 (u fj u { j) 

V ij * s 1j /(8 kl s kl )1/2 > ' ,/2 




1/2 


• u ij m 9 ij^°ki°ki 


1j a 1j / ' l0 krkl J 
<»n/(®kA 


, 1/2 


with Z (o) * Z 0 ; W p » ; W p (o) « 0 ; a( 0 ) » 0 


b. Oirectional Hardening 

®1J ’ m 2 (Z 3 u 1j * 9 1j ) *p * Vl 


tSkl Ski ) 
2l 


1/2 


•a 


with Z D « ^1 J U 1 j ; 2 ° (o) ■ 0 , 8^(0) • 0 


Material Constants: D 0 , Z 0 , Z] , Z 2 , Z3, , m2, at] 

At, A 2 * r lt r 2 , n, and elastic constants 


> - 


In most cases can set: = r^, A^ * A^, Z Q = Z ,, 





TABLE 3.2 


MATERIAL CONSTANTS IN BODNER-PARTOM CONSTITUTIVE MODEL 


Constants 


h 

h 

Z 3 

m i 

m2 

n 

A: 

h 

r i 

r 2 

°1 


Descriptions 

Limiting shear strain rate, see'* 

Initial value of the isotropic 
hardening variable, MPa 

Limiting (maximum) value of Z*, MPa 

Fully recovered (minimum) value of Z*, MPa 

Limiting (maximum) value of Z^, MPa 

Hardening rate coefficient of Z*, MPa'* 

Hardening rate coefficient of Z°, MPa'* 

Kinetic parameter 

Recovery coefficient for Z*, sec"* 

Recovery coefficient for Z* 1 , sec"* 

Recovery exponent for Z* 

Recovery exponent for Z° 

Limiting coefficient for non-proportional 
hardening effect 
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represent material behavior over a wide range of strain rates including the 
very high rates, i.e., about 10 5 sec - *. 

The treatment of directional hardening by Bodner and his associates 
has been based on the general anisotropic form of the flow law rather than on 
the "back stress" concept. An incremental isotropic form of the flow law has 
been proposed as an approximation which would simplify the numerical computa- 
tions. In the incremental isotropic equatijn, the scalar coefficient is a 
function of both the isotropic and directional hardening variables which de- 
pend on the loading history. The isotropic hardening variable controls the 
cyclic hardening or softening characteristic. 

The scalar measure for hardening is taken to be the plastic work rate. 
However, the overall theory is not dependent on this point and other measures, 
e.g., plastic strain rate or a function of the plastic strain rate history, 
are admissible within the context of the theory. For directional hardening, 
the stress is used as the directional index of hardening to avoid cross soft- 
ening effects. The plastic strain rate had been used previously for that pur- 
pose and most other theories still rely on that variable. 

Multiaxial, non-proportional hardening is taken into account by intro- 
ducing a load history dependent variable, a, into the evolution equation for 
isotropic hardening. As can be seen from the expression in Table 3.1, the 
parameter a increases with non-proportional loading history from zero to a 
limiting value of a^. The saturation value of Z* under such loading is then 
(Zj+ aj Z3). Non-proportionality of loading path is measured in terms of an 
angular parameter e defined by either of the following expressions: 


! 


.1 


f 

1 


e = cos 
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(3.1) 


V 



0 


t 


t \ 


or 0 - cos -1 (U^. U^) (,.2) 

where U^j, U^., V^j and V.^., as given in Table 3.1, are the directional 
cosines of the current values of a^, s^j, o^-, and 6^, respectively. 

3.2 Walker Model 

Two internal variables are used in the Walker theory. The equilibrium 
stress fljj represents nonlinear kinematic hardening and accounts for the 
Bauschinger effect, while the drag stress K represents isotropic hardening and 
accounts for cyclic hardening of the material. 

The growth law .or the equilibrium stress, n. it is the sum of two com- 
ponents, a^j Q and n^. a.j 0 is the initial equilibrium stress offset which 
accounts for the difference in the tensile and compressive yield strength. 
The growth law for the second component, n.jp contains hardening, dynamic 
recovery and thermal recovery terms. At high strain rates thermal recovery is 
small in comparison to the dynamic recovery and the equilibrium stress 
saturates to a limiting value at large plastic strain that is independent of 
the strain rate. The drag stress, K, is assumed to be a function of only the 
accumulated effective plastic strain, R.* The isotropic hardness is initially 
- Kp an d saturates at a value of with increasing plastic strain. Cyclic 
hardening or softening can be modelled using this feature. The kinetic rela- 
tionship that relates the nonlinear strain rate to the excess stress, S^j- 
was assumed to be a power law expression. This form of the model [3] has 
been found to capture most of the behaviors noted in testing of the solid 
solution strengthened alloys such as Hastelloy-X. 


R = 


ft 
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The Walker model, as described in [3] and summarized in [1], has been 
significantly modified for describing the constitutive behavior of BlOOO+Hf. 
Modifications to the unified constitutive model to capture some of the fea- 
tures of the B1900+Hf deformation behavior are indicated: 

(1) Form of the Kinetic Relationship - At low to intermediate 
temperatures the limit stress in a monotonic tensile test 
does not vary with strain rate to a power but exponential- 
ly. A more appropriate equation relating the inelastic 
strain rate and the effective stress is the exponential 
function: 


•P _ i ovn- 

>i 

[2/1 s 
'3 1 2 H 


e - e ex P| 



K _i 


(2) Form of the Drag Stress Relationship, K - Cyclic hardening 
or softening observed during the uniaxial testing is simu- 
lated by assuming that the isotropic hardening variable, 
K, varies with accumulated inelastic strain, R. Addition 
of a similar term, K 3 x exp (- ngL), can approximate the 
cyclic hardening observed under nonproportional loading 
cycles. L, a measure of the nonproportionality of the 
strain path, is defined as the difference between the ef- 
fective strain rate and the rate of effective strain, as 
shown in Eq. 3 of Table 3.3. ng and K 3 are shape and 
level parameters, respectively, and can be determined from 
a single nonproportional multiaxial loading test at each 
test temDerature. 






V 




\<a. \ 


^S: 




(3) Form of the Equilibrium Stress Relationship - The equili- 
brium stress is now assumed to be the sum of three terms, 
fl i jl + ^ijo* an< * an a ^ 1 ^^ ona ^ term, introduced to 

predict linear hardening at high strain rates. The ex- 
pression for 


• j2 


n u e 


P 

ij 


‘ n ij2 (n < 


R + 


10 ' 


(3.4) 


For plastic strain rates sufficient to saturate a. . + 

I Jv 

but insufficient to saturate 8.^, varies 
linearly ^th the plastic strain rate. Assuming a con- 
stant strain rate tensile test, the flow stress will 
therefore increase directly with a.j 2 in a linear fashion 
for small strains. With increasing strain and for high 
strain rates, will eventually saturate to a value 
n^/ng. With increasing strain and at low strain rates, 
the rate of strain hardening in can be minimized by a 
judicious selection of hjq, representing thermal recovery. 
A modification to the evolutionary expression for is 
also required to capture the inverse strain rate sensiti- 
vity observed in the B1900+Hf monotonic tensile data at 
760°C. The dynamic recovery term is now assumed to be a 
function of the strain rate, vis: 


(n 3 + n 4 exp ( -n 5 | log (-j-) \ ) 


(3.5) 
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By judicious selection of the constants n 4 and n 5 the flow 
strength can be made to decrease either side of the strain 
rate Rq. 

(4) Thermal Effects - The constants in Walker's unified theory 
are generally functions of temperatures, T. If a material 
is subjected to a strain and temperature change simultan- 
eously, the incremental change in a state variable, SV, 
will be: 


dSV aSV aSV dT 

dt at aT dt 


(3.6) 


The Walker theory has been modified to include the rate of change of 
temperature terms. 

Table 3.3 presents the modified Walker theory employed by PWA to cor- 
relate the B1900+Hf data. The modified Walker model contains 15 material con- 
stants which are all temperature dependent. Generally, fewer will be needed 
(about 10 for B1900+Hf) in an actual material representation. 
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4.0 UNIAXIAL EVALUATION OF CONSTITUTIVE MODELS 


4.1 Experimental Procedures 

Five different mechanical testing modes were used for developing the 
uniaxial constitutive data of B1900+Hf in the present program. They included 
tensile, creep, isothermal cyclic, stress relaxation, and thermomechanical 
fatigue (TMF) tests. A large quantity of constitutive information was also 
generated as part of the data acquisition efforts in the HOST Isotropic 
Fatigue Contract (NAS3-23288). Results of these two programs were integrated 
to form a large data base for constitutive model evaluation. 

All the isothermal, uniaxial tests in the present program were con- 
ducted in a closed-loop, servo-hydraulic test machine at SwRI. Figure 4.1 
shows the configuration of the test specimen and the experimental set-up. 
Strain was measured using an externally mounted extensometer with ceramic 
reach-rods attached to the specimen. Induction heating was used for providing 
the elevated temperature; temperature variations over the specimen gage sec- 
tion are within ASTM specifications for short-term tests. The specimen de- 
signs utilized in various constitutive tests are schematically illustrated in 
Figure 4.2. Figure 4.3 shows photographs of: (A) an isothermal constitutive 
specimen, (3) a thermomechanical constitutive specimen, and (C) a biaxial spe- 
cimen. 

Tensile tests were performed at various strain rates ranging from 
4xl0' 2 to 4xl0“ 7 sec -1 at 25, 649, 760, 871, 982, and 1093°C. A summary of 
the tensile test matrix can be found in last year annual report [1!. Creep 
tests were conducted at the latter five temperatures. 
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FIGURE 4.2. SPECIMEN DESIGNS UTILIZED IN VARIOUS CONSTITUTIVE TESTS. 
(All units are in inches.) 
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FIGURE 4.3. PHOTOGRAPHS OF (A) ISOTHERMAL CONSTITUTIVE 
SPECIMEN, (B) THERMOMECHANICAL CONSTITUTIVE 
SPECIMEN, AND (C) BIAXIAL SPECIMEN. 


TABLE 4.1 


■ ^ Hat. a . 

u f* 


ISOTHERMAL CYCLIC CONSTITUTIVE TESTS 


M 


i 

r 

i 

! 

| 



Strain Range 



Temp(°C) 

R 

e(sec _1 ) 

.25% 

.4% 

.6J5 

.8% 

0 

Effect of Temperature and Strain Rate 





427 

0 

4xl0’ 3 

X 

X 

X 

X 


538 

0 

4xl0“ 4 

X 

X 

X 

X 


649 

0 

4xl0" 4 

X 

X 

X 

X 


760 

-1 

4xl0 -3 

X 

X 

X 

X 


871 

-1 

X 

o 

1 

4* 

X 

X 

X 

X 


982 

-1 

4xl0 -4 

X 

X 

X 

X 


1093 

-1 

4xl0' 5 

X 

X 

X 

X 

0 

Effect of Mean Strain 







538 

-1 

4xl0' 4 

X 

X 

X 

X 


760 

0 

4xl0~ 3 

X 

X 

X 

X 


760 

— CD 

4xl0“ 3 

X 

X 

X 

X 


871 

— CD 

4xl0 -4 

X 

X 

X 

X 


1093 

— OO 

4xl0~ 5 

X 

X 

X 

X 

0~ 

Strain Hold to Probe for Equilibrium Stress 





760 

-1 

4xl0~ 4 

X 

X 

X 

X 


871 

-1 

4xl0“ 4 

X 

X 

X 

X 


1093 

-1 

4xl0" 4 

X 

X 

X 

X 

r 

Effect of Stress Hold 







Temp(°C) 

R 

a (MPa/sec) 

taoj 

tAOp 

tAo^ 

± Aa 4 





(MPa) 




760 

-1/T 

59.3 

t 350 

± 550 

± 650 

± 750 


871 

-1/TC 

55.2 

± 310 

± 400 

i 500 

± 600 


982 

-1/C 

51.0 

i 241 

± 345 

— 

— 


T - Stress hold at peak tension for 120 sec. 

C - Stress hold at peak compression for 120 sec. 
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Isothermal cyclic tests were conducted either with or without dwell. 
As indicated in Table 4.1, the cyclic constitutive tests were performed at 
four different strain range levels (±.25%, ±.4%, t.6% and ±.8%). Testing was 
performed at R ratios (R = ^p/emax) of °« -1. and -■=■> at temperatures ranging 
from 427 to 1093°C and at strain rates of ±4x10“^ to ±4:.10"^ sec“^. Cyclic 
tests with dwell were conducted at a loading rate of approximately 58.6 MPa/s. 
The load was held constant for two minutes when the prescribed maximum or 
minimum stress was reached. For determining equilibrium stress values, stress 
relaxation at various strain levels of the first cycle hysteresis loop and of 
the saturated one was measured. The time of strain hold was 30 seconds to 2 
minutes. The stress relaxation test conditions are also shown in Table 4.1. 

The thermomechanical fatigue (TMF) loading histories were run with a 
sinusoidally varying mechanical strain and temperature with a cycle period of 
sixty seconds (Figure 4.4). The strain for three of the tests varied from a 
nominal -0.4 percent to +0.4 percent, and the temperature varied sinusoidally 
from 538° to 982°C. For the first TMF loading eyrie, the temperature and 
strain varied in phase, i.e., the maximum temperature occurred at the same 
time as the maximum strain (Figure 4.4a). For the second loading case, the 
temperature and strain loads were applied 180° out of phase, and in this case 
the maximum temperature coincided with the minimum strain (Figure 4.4b). The 
third TMF cycle was a repeat of the second case but with the addition of a 
sixty second constant strain hold at -0.4 percent strain. 

4.2 Evaluation of Model Constants 

4.2.1 Bodner-Partom Model 


A systematic procedure for determining the material constants in the 
Bodner-Partom model is outlined in Figure 4.5. The procedure is developed by 


Tensile Data 










Tensile Data and/or 
Creeo Oata 



a -o n 
c o 


FIGURE 4.5. A FLOW CHART SHOWING A SYSTEMATIC PROCEDURE 
FOR EVALUATING THE MATERIAL CONSTANTS IN 
THE BODNER-PARTOM THEORY. 
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recasting the Bodner-Partom equations Into convenient forms such that the 
material constants can be evaluated from the experimental data by separate 
consideration of basic physical effects such as hardening, strain rate depen- 
dence, and thermal recovery. The theoretical basis of the method is outlined 
here. Excluding the elastic constants, a total of 13 material constants must 
be evaluated from pertinent experimental data. In most cases, including non- 
proportional loading paths, only nine model constants need to be determined. 
Most of the constants, tabulated and briefly discussed in Table 3.2, can be 
evaluated from monotonic tensile data. Creep data would also be required if 
slow strain rate tensile data were not available. 

As shown in Figure 4.5, step (1) in determining the material constants 
in the Bodner-Partom model is to develop a set of work hardening rate data 
from the experimental stress-strain results. A simple computer program has 
been developed for this purpose. Using the monotonic tensile stress-strain 
data and the elastic modulus as inputs, the computer program calculates the 
plastic strain component by subtracting the elastic component from the total 
strain. The stress is then expressed as a polynominal function of the plastic 
strain. 


2 m 

o = a Q + a l£ p + a 2 (e p ) + + a m ( £ p ) (4.1) 

The coefficients a Q , aj, a 2 , ... a m are computed using the least-square 
method. The work hardening rate, y = do/dW^ = (l/o) (do/d £ p ) , can then be 

calculated as: 


where a is obtained from Equation 4.1. The number of terms required In Equa- 
tion 4.1 is usually obtained by trial and error. For B1900+Hf, the cubic 
polyncminal (m - 3) is generally sufficient, but a fourth-order (m = 4) poly- 
nominal has also been used occasionally. 

Step (2) is to evaluate the constants nij and from the plot of y 
versus o (Figure 4.6). The Bodner-Partom model can be written for uniaxial 
stress case as: 

a = L 1 Z (W p ) (4.3) 

n *l/2n 

with Lj = (2 In (2D Q /v / 3 e P ) ] 

Z * Z 1 + Z D 

where Z 1 and Z° are the Isotropic and the directional hardening variables. 

In the absence of thermal recovery and on the physical basis that Z D 
saturates much more rapidly than Z* for most materials, i.e., m? > m^, the 
following approximations to Equation (4.3) can be made in the region of small 
plastic strains. 


Z 1 = (o/Lj) - Z° * (o/Lj) - Z 3 
Z D = (o/Lj) - Z 1 . (o/L 1 ) - Z Q 

where the additional terms are defined in Table 3.2. These approximations are 
then substituted in the general expression for y = do/dWp obtained from 
Equation (4.3) and the evolution equations for Z 1 and Z D , namely 
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FIGURE 4.6. A SCHEMATIC SHOWING THE ISOTROPIC AND DIRECTIONAL HARDENING 
COMPONENTS IN A y-o PLOT. 



Y 


(4.4) 


- 4 L 1 (Z 1 -Z I ) * m 2 (Z 3 -Z°)j 

Upon reduction, the expression for y becomes 

r = L 1 [m 1 (Z 1 +Z 3 ) + m-(Z o +Z 3 )J - (mj+m^o (4.5) 

This indicates that a plot of y against a should be linear with slope (m^ + 
m^) in the range of small plastic strains. 

At the larger strains, Z° fully saturates to the limiting value Z 3 so 
that setting Z° = Z 3 and again using the above approximation for Z 1 gives 

y = L^m^(Z^+Zj) - ay (4.6) 

which is a linear relation with slope ny These relations indicate that the 
y-o function should be essentially linear at both the low and high plastic 
strain regions. 

Plots of y -a in Figure 4.7 indicates the bilinear behavior. For tem- 
peratures at and below 871°C, the effects of the strain rate and temperature 
are to shift the y-a curves along the horizontal axis without causing any sig- 
nificant changes in the slopes of the curves. This hardening behavior can be 
rationalized on the basis that both the isotropic and the directional harden- 
ing components are present and they are not significantly affected by static 
tnermal recovery at these temperatures. On this basis, the constants mj and 
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Stress, <r, MPn 


FIGURE 4. 


(b) 


. THE y-o CURVES OF 31900+Hf (A) AT 871 °C SHOWING 
SHIFTS TO HIGHER STRESSES AS THE STRAIN RATES 
ARE INCREASED, AND (8) AT A GIVEN STRAIN RATE 
SHOWING TEMPERATURE SHIFTS ALONG THE a-AXIS AND 
SLOPE CHANGES FOR TEMPERATURES ABOVE 982°C. 











m 2 can be evaluated from a y -a plot. As shown in Figure 4.6, the initial 
linear slope of the y-o curve for B1900+Hf at 871°C is given by the sum of m^ 
and m 2 , while at high stress values the slope is then given by m^. The extra- 
polated intercept at the a axis for y = 0 is o $ (the saturation stress) which 
is given by: 


o $ = L^(Z^+z_ 3 ) (4.7) 

The measured values of m^, + m 2 , o , and the calculated values of m 2 for a 

number of specimens tested at various temperatures are summarized in Table 
4.2. 

Step (3) is to evaluate the constants n and D Q in the kinetic equa- 
tion. D 0 is the limiting inelastic strain rate in shear and is usually taken 
to be lxlO 4 sec -1 unless high strain rates are being considered. The constant 
n is evaluated from the stress dependence on strain rate at a constant value 
of Z. Equations 4.3 and 4.4 are rearranged to give the following expressions: 

h no ^ In S + In Z (A. 8) 

and -pj- = - ; Z = constant; W p = constant (4.9) 


where S = 2 Zrc (20 o //3 e P ) = Lj' 2n 


(4.10) 


-i v-e 


*1 

\ 
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EVALUATIONS OF BODNER-PARTOM MATERIAL CONSTANTS FOR B1900+Hf 
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Thus, n can be evaluated either from the slope of a plot of Ina versus InS at 
a constant value of Z (Wp) or from Equation (4.9). When static thermal recov- 
ery is not important, the stresses for evaluating the parameter n could be 
( 1 ) the saturation stress o $ at different strain rates where a $ corresponds to 
the limiting value of Z or (Z^ + Z 3 ) which can be obtained by extrapolating 
the y -0 curve to the stress at which y is equal to zero (Figure 4.6); (2) the 
peak stresses of saturated cyclic hysteresis loops at different strain rates 
and at a constant value of large strain amplitude for which both Z* and Z^ 
saturate; and (3) the flow stresses at different strain rates and at constant 
values of plastic work. The n values for B1900+Hf at 871°C and 539°C have 
been calculated based on method (1) using Equation (4.9). The results are 
shown in Table 4.2. 

It is important to note that stress values under steady state condi- 
tions during creep or tensile straining do not always correspond to the same 
value of Z. In addition, plastic work does not have a one-to-one relationship 
with Z at elevated temperatures for which static thermal recovery is signifi- 
cant. Consequently, stresses under steady state conditions involving thermal 
recovery should not be used in conjunction with Equations 4.8 and 4.9 for 
evaluating n. Once ! ne n value is evaluated for a lower temperature at which 
static thermal recovery is insignificant, the constants Z^ and Z 3 can be com- 
puted (see Step 4). Since both Zj and Z 3 are considered to be temperature- 
independent, the n value at high temperatures can be obtained by fitting a 
stress-strain curve at a strain rate above 10"^ sec"* or by using Equation 
(4.7) if a $ is known. The n values for B1900+Hf at 982 and 1093°C have been 
obtained using the former procedure. 


Step (4) is to calculate the sum of and Z 3 . Using the saturation 
stresses obtained in Step (2) ai*d the value obtained in Step (3), the sum 
of Zj and Z 3 is calculated using Equations (4.7). The results are shown in 
Table 4.2. 

Step (5) is to evaluate the relative contributions of Zj and Z 3 , the 
limiting values of the isotropic hardening, and the directional hardening 
variables, respectively, in the total hardening. From Equation 4.6 and Figure 
4.6, it can be shown that is given by: 

m l4 (Z l + Z 3* + m 2 L l (Z o + Z 3> „ 

°b = STTin < 4 - u > 


Equation 4.11 can be rearranged to yield 

, _ ( m i +m 2 ^ °b ‘ m l°s " m 2 °o 
L 2 m 2 L 1 


(4.12) 


with 


4 z o (T) 


(4.13) 


where o Q is taken as the . 2 % offset yield stress, i.e., an initial value of 

Z. The values of o. and 0 and the calculated values of are summarized in 

bo j 

Table 4.2. Also, it is assumed Z£ = Z Q for most materials (Step 6 ). For spe- 
cimens subjected to full thermal recovery, Z Q is best estimated using slow 
rate tensile results. Once Z 3 and Z^ are obtained, Z 0 can be refined by fit- 
ting the model calculation to the tensile stress-strain curve (Step 6 ), and Z£ 
is again set to equal Z Q . 

The effect of static thermal recovery on the work hardening rate can 
be included in a similar manner. Assuming A^ = A 2 = A and r^ = r 2 = r, the 


expression for y including the recovery term can be shown to have the form 
(see Appendix): 


Y 


— L ^ [fiij (Z^+Z^) n^fZg+Z^) 1 - (m^+m^) 0 


A(L 1 Z 1 ) 


1-r 


oe 



(4.14) 


Equation 4.14 is reduced to Equation 4,6 when A is zero; this occurs when sta- 
tic thermal recovery is unimportant. Equation 4.14 also indicates that the 
dependence of r on a is altered when static thermal recovery is present. 
Thus, it is possible to identify the temperature regime at which static ther- 
mal recovery is important by comparing y-a data over a range of temperatures 
(e.g., see Figure 4.7). 

The final step (Step 7) is to evaluate the recovery constants. For 
B1900+Hf, rj is assumed to be equal to r 2 (r^ = r 2 = r) and A^ is also assumed 
to be equal to A 2 (Aj = A 2 = A) for the following reasons: (1) there is a 
need to reduce the number of constants and to simplify the procedures for 
evaluating these constants, (2) the r-o curves in Figure 4.7 suggest that they 
need not be different an^‘ might be difficult to distinguish if they were, and 
(3) no physical mechanisms suggest that they should be different. The recov- 
ery constants can be evaluated from either the steady state stress (o c ) and 
associated strain rate (i c ) of a slow rate tensile test or the steady state 

creep rate (e ) and corresponding stress (a ) of a creep test. In both cases, 
c c 

the stresses and the plastic strain rates do not change with time or deforma- 
tion. Under these conditions, the work hardening parameter, y» is zero and 
Equation 28 can be rewritten as: 



0 


(4.15) 


q - A(L 1 Z 1 ) 1 - r (o c -a Q ) r = 


where q = + n» 2 (Z 0 + Z 3 )l - (m 1 +m 2 )a c } 


The recovery exponent r can be obtained from Equation (4.15) as: 


d log q 
d log (<3^-a 0 ) 



(4.16) 


Figure 4.8 shows data of q vs ° c ~° 0 in a log-log plot for B1900+Hf at 871, 
982,, and 1093°C. The results Indicate that the values of r are equal to 2 at 
these temperatures. Once r Is obtained, the recovery constant A can be calcu- 
lated using Equation (4.15). Alternately, A can also be evaluated by consid- 
ering individually the conditions for which the isotropic hardening and the 
directional hardening rates are zero. In general, they both give the same re- 
sult. If different A values are obtained, the one which gives the best fit to 
the tensile stress-strain curve at e = 4x10"^ sec”* is used. 

A summary of the material constants of B1900+Hf using this set of pro- 
cedures is shown in Table 4.3 as well as in Figure 4.9. There are only four 
temperature-dependent materials constants in the Bodner-Partom model, namely, 
n, Z Q , Z 2 and A. Since Z 2 is generally taken to equal Z Q , only Z Q , n and A 
need to be evaluated as a function of temperature. For the B1900+Hf alloy, 
Z Q , n, and A are temperature-dependent only at temperatures above 760°C. The 
dependence of l Q or temperature can be described by the following two empiri- 
cal expressions: 
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FIGURE 4.8. THE SLOPES OF THE q VS. o c -o 0 CURVES FOR 
B1 900+Hf ARE APPROXIMATELY TWO AT 871, 
982, AND 1093°C. 



Temperature- Independent Constants 


hmA 



8.650 x 10 - 17.58 T + 2.321 x 10 T - 3.464 x 10 T MPa with T in 









Ni 



(a) 



.(b) 


FIGURE 4.10. THE ELASTIC AND SHEAR MODULI OF B1 900+Hf 
AS FUNCTION OF TEMPERATURE: (A) ELASTIC 

MODULUS, (B) SHEAR MODULUS. 







2700 MPa, T < 838°K 


(4.17) 


Z 


o 


and Z Q = - 4950 + MPa, T > 863.4 ft 

The empirical relationships relating n and temperature are: 


n = 1.055 for T < 863. 4°K (4.18) 

and n = -1.0524 + — 4 — 7 - 4 for T > 863. 4*K 

The constant A in the thermal recovery term is best expressed as an exponen- 
tial function of temperature. The choice of the functional form is suggested 
by Friedel's recovery model [ 10 1 and based on the consideration that thermal 
recovery is controlled by diffusion and should have a similar functional form. 
As shown in Figure 4.10, the constant A for 81900+Hf can be expressed as: 


A = 1.319xl0 9 exp J sec 1 


(4.19) 


where T is in °K. It should be noted that at temperatures below 760°C, the 
value of A obtained from Equation (4.19) is negligibly small so that it is 
essentially zero and thermal recovery is unimportant. The elastic and the 
shear moduli obtained from the slope of the tensile and the shear strain- 
strain curves were also expressed as a function of temperature by fitting data 
to a third order polynominal equation. The expressions for the elastic and 
the shear moduli of B1900+Hf, shown in Table 4.3, are compared with the ex- 
perimental data in Figure 4.10(a) and (b), respectively. 

4.2.2 Modified Walker Model 

Two strategies have been developed for determining constants in the 
Walker constitutive theory. The first approach requires simple tensile or 
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creep testing supplemented with limited cyclic data to establish, in a 
straightforward manner, the variation with temperature of the model constants. 
Test and an?' sis requirements are minimized but there is nc guarantee, a 
priori, that the resultant model will reproduce material behavior under com- 
ponent relevant strain cycles. The second approach, which is described in 
[1), proceeds from the notion that the most appropriate test for determining 
constants is that which most closely represents the intended application. 
This approach necessarily involves considerable trial and error before test 
data and predictions match. The use of optimization techniques has been found 
to be of some value in this regard. 

All the Walker model constancs for B1900+Hf in this report have been 
determined using the first approach. An outline of the systematic procedure 
to establish the Walker model constants is presented in Figure 4.11. The data 
requirement at each temperature consists of a series of monotonic tensile 
tests at different strain rates and/or creep tests at different stress levels 
and fully reversed cyclic tests at one, but preferably two strain rates. The 
monotonic tensile stress-strain results at high strain rates allow n^» n g , 
and n^Q to be estimated directly, is seen to be approximately equal to 

the initial value of the strain hardening rate, do/de p . The difference 
between the limit and the yield strength determines the ratio, n^/ng. 

n 10 fi 21iml is selected that ( and therefore n 211m^ 1s sma11 except for 
the strain rates e* 'u i linear hardening is noted. Iterative analysis is 
required to estabii-.i the remaining constants. A plot of the limiting or 
large strain value of the excess stress ( S^j - - ^lim^ against the 

log of Ip should be linear with a slope of and one cycle intercept beta, 
but ^2lim’ as we have seen » and n liim are typically functions of strain rate. 
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FIGURE 4.11. A SCHEMATIC PROCEDURE FOR EVALUATION OF THE 
MATERIAL CONSTANTS IN THE EXPONENTIAL FORM 
OF THE WALKER THEORY. 
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Determination of the constants through ng begins by plotting the quantity 
(°1 im~ n 21 in/*p^ versus the log of Cp. A straight line indicates that strain 
aging (n^ and ng) and thermal recovery (ng) effects are small and is a 

constant, independent of the plastic strain rate. A sigmoidal curve indicates 
that thermal recovery (ng) is present and a local peak in the o^ m - 
(ep) versus log of £p plot at a strain rate of £p Q indicates that strain aging 
effects must also be modelled. At limit conditions, = 0, and 


n llim 


n 3 /,T1 2 + n 4/ n 2 e 


— Hg | lo g e p /C pQ l 


+ "eVp 


At high strain rates, significantly removed from £p Q , the second and 
third terms in the denominator can be ignored and the ratio n^/n^ can be set 
equal to a provisional value of Kj can be determined as the slope of 

(°lim " n lliml ~ n 21iml^ versus lo 9 of in this region of the curve. At 
very low strain rates, the third term in the denominator dominates and the 
ratio ng/r^ can be determined from the relationship 


n g/ n 2 = 


£ P4 

s 11im4 


where a^.^ is the saturated back stress at a plastic strain rate e p4 

required to linearize the excess stress (o,. - a 01 . - a... ) versus log of 

Ep plot. Finally, n 5 and the ratio 114/112 are selected to match the variation 

of with due to strain aging effects in the vicinity of e . n^is now 

i nm p po2 

determined to provide an optimum fit of the monotonic tensile curve (variation 
of c p with 0 - n^p)) at high strain rate. The constants n^» n 4 and ig ars, 
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FIGURE 4.12. WALKER MODEL CONSTANTS FOR B1900+Hf AS FUNCTION 
OF TEMPERATURE. 





therefore, determined for the provisional value of 8 can be solved 

directly, completing the trial determination of the constants. 

The material constants for the modified Walker model for 81900+Hf are 
presented in Table 4.4 as well as in Figure 4.12 as a function of temperature. 
Isothermal tensile and a limited number of one way cyclic test results were 
used to determine these values. Creep, dwell fatigue, and TMF data were not 
considered in the selection of the model constants. 

4.3 Correlation of Theoretical and Experimental Results 

4.3.1 Uniaxial Tensile Behavior 

Figure 4.13 shows the monotonic tensile results and corresponding 

model predictions at three temperatures and a fixed strain rate of 8.3 x 10"^ 

sec'^. Both models provide a good fit of the data. The effect of strain rate 

on the 0.2% yield stress at 760, 871, 982, and 1093°C is shown in Figure 

4.14. The Bodner-Partom model predicts an increase of yield stress with 

strain rate at all temperatures. The thermal recovery term is absent in the 

760°C calculations but is dominant at 1093°C. The change in the slopes of the 

Bodner-Partom model calculations at 871 and 982°C represents a transition from 

-5 -1 

a thermal recovery dominated regime at slow strain rates (e < 1x10 sec ) 
to a hardening dominated regime with relatively little thermal recovery at 
high strain rates (e > lxlO'^sec'^) . By including a strain-aging term in the 
equilibrium stress, the Walker model is able to correlate the decrease in 
yield strength with strain rate that is observed at 760°C at e > lxlO'^sec"^. 

The temperature depender . of the .2% e P yield stress at various 
strain rates is shown in Figure 4.15. The Bodner-Partom model reproduces the 
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FIGURE 4.13. COMPARISON OF EXPERIMENTAL AND CALCULATED STRESS 
STRAIN CURVES OF B1900+Hf AT THREE TEMPERATURES. 
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FIGURE 4.14. THE EXPERIMENTAL AND CALCULATED .2% e 1 " YIELD STRESS OF B1900+Hf 
AT VARIOUS TEMPERATURES AS FUNCTION OF STRAIN RATE. 









experimental data rather well. The expp. imental and calculated yield stresses 
are seen to be relatively constant and independent of temperature at T < 760°C 
but decrease rapidly with temperature at T > 760°C. 

4.3.2 Creep Behavior 

The variation of the steady-state creep rate with test temperature and 
stress level is presented in Figure 4.16. Note that predictions obtained 
using the Bodner-Partom model give somewhat better agreement with the experi- 
mental data than the Walker model at 760°C and at low strain rates at 871°C. 
This finding is somewhat surprising considering the good fit of the tensile 
data at these temperatures obtained using the Walker model. 

Figure 4.17 shows creep strain vs time at various stress levels at 
871°C and 1093°C. The unified models in the present form predict primary i d 
steady state creep regions but not a tertiary region. Comparison of model 
calculations and experimental data indicates that the creep strain at a given 
time is predicted fairly accurately by the Bodner-Partom model. It should be 
noted that the creep curves generally show substantial scatter so that eval- 
uation of model accuracy on the creep response should take this into account. 

4.3.3 Uniaxial Cyclic Constitutive Behavior 

Representative predicted and experimentally observed saturated hyster- 
esis loops for fully reversed strain cycling at 538°C, 871°C, and 1093°C are 
shown in Figure 4.18, 4.19, and 4.20, respectively. The shape of the hystere- 
sis loops are predicted fairly well by both models at 538 and 871°C. At 
1093°C, both models predict hysteresis loops which tend to be more square than 
the experimental ones. Furthermore, the Walker model predicted hysteresis 
loop also shows bi-1 inear regions. 
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The cyclic stress-strain curves were constructed from both the experi- 
mental and theoretical results by plotting half of the stress ranges of the 
*u11y-reversed saturated hysteresis loop vs half of the imposed strain range. 
Comparisons of model predictions and experimental data for four different tem- 
peratures (760, 871, 982, and i093°C) arc summarized in Figure 4.21. Note 
that the Bodner-Partom model overpredicted the stress ranges at 871 and 982°C. 
However, it is worthy to note that none of the cyclic data were used for 
determining Bodner-Partom constants. 

The effects of imposed compressive and tensile mean strains on the 
cyclic constitutive behavior of B1900+Hf are summarized in Figures 4.22, 4.23, 
4.24, and 4.25 with corresponding predictions using both theories at 1093°C, 
871°C, TSO’C, and 538°C, respectively. For purposes of comparison the half 
stress range, half cyclic plastic strain range, and mean cyclic stress at both 
the first and sixth cycle are plotted versus the half strain range for R 
ratios of 0, -1, and -=>. Two important observations in Figures 4.22-4.25 are: 
(1) the experimental and theoretical saturated cyclic stress-strain curves 
( ao/ 2 - Ae/2 curves) appear to be unique at a particular strain rate and tem- 
perature and are independent of the R ratios, and (2) both unif,ad models pre- 
dict a drift in the mean cyclic stress which is not always observed in the ex- 
perimental data. 

At 1093°C, Figure 4.22, both models yield adequate fits of the stress 
range and cyclic plastic strain range, but the Bodner-Partom model predicts 
the r.tress amplitude more accurately in the lower strain range. Both consti- 
tutive models match the sixth cycle response of B1900+Hf at 871 C quite well 
in all three aspects (Figure 4.23). The predicted first cycle mean stress for 
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FIGURE 4.21 


COMPARISON OF EXPERIMENTAL AND PREDICTED CYCLIC 
STRESS-STRAIN DATA OF BI900+Hf OBTAINED BY PLOT 
TING HALF OF THE STRESS RANGE VS. HALF OF THE 
STRAIN RANGE OF THE STABLE HYSTERESIS LOOP. 
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FIGURE 4.23. PREDICTED AND EXPERIMENTAL CYCLIC STRESS- STRAIN 
DATA, MEAN STRESS, AND PLASTIC STRAIN RANGE OF 
B1900+Hf UNDER TWO MEAN STRAIN LEVELS AT 871°C. 
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FIGURE 4.25. 


PREDICTED AND MEASURED CYCLIC STRESS-STRAIN 
DATA, MEAN STRESS, AND PLASTIC STRAIN RANGE 
OF B1 900+Hf AT NO OR A COMPRESSIVE MEAN 
STRAIN AT 538°C. 




compression cycling is overpredicted usina the Walker model, however. At 
760°C the Bodner-Partom model prediction agrees with the data experimental 
quite well in all respects (Figure 4.24). Walker model predicts essentially 
the same stress range but slightly greater cyclic plasticity than the Bodner- 
Partom model and a pronounced drift in the cyclic mean stress that was not 
observed in the data. First cycle test data at 538°C, Figure 4.25, is in good 
agreement with the predictions obtained using the Walker model but, as noted 
above, the mean cyclic stress for eccentric strain cycles is predicted to 
drift appreciably. 

Figure 4.22-4.25 also indicates that the Bodner-Partom model shows re- 
latively small mean stress drift when the plastic strain range is small, but 
becomes increasingly large when the plastic strain range becomes relatively 
large. The mean stress drift as a function of cycle number of both the 
Bodner-Partom and Walker model is compared with experimental data at rela- 
tively large plastic strains in Figure 4.26. The results indicate that both 
models predict mean stress relaxation in less than ten cycles, whereas the 
mean stresses are maintained in the experiment during the whole test period. 
It should be noted, however, that the mean cyclic stress in these strain 
cycles is small compared to the yield stress under the test conditions. 

4.3.4 Stress Relaxation at Strain Hold 

Strain holds were performed at various strain levels of hysteresis 
loops for 2 minutes to probe for the equilibrium stresses associated with cyc- 
lic deformation of B1900+Hf. Figure 4.27, which presents typical results of 
the equilibrium stress probe test at 760°C, indicates that stress relaxation 
occurs during strain hold but subsequent rapid hardening increases the stress 
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FIGURE 4.26. DRIFT OF MEAN STRESS AS FUNCTION OF CYCLE NUMBER. 
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FIGURE 4.27. STRESS RELAXATION DURING STRAIN HOLD IN THE EQUILIBRIUM 
STRESS PROBE TEST. 





back to the previous value. Cyclic hardening is evidenced since the stress 
range of the first loop is less than that of the saturated loop which is the 
22nd in this case. The stress ranges of the saturated loops with strain hold 
are compared with loops without strain hold in Figures 4.23 and 4.24. The 
comparison reveals no significant difference, indicating the stress ranges of 
the saturated loops with strain hold are not affected by the prior deformation 
history. 

Stress relaxation data associated with holding at the peak strain of 
various saturated hysteresis loops are shown in Figures 4.28-4.29, together 
with the Bodner-Partom model predictions. In general, the agreement between 
theory and experiment is good for the three test temperatures which include 
871, 982, and 1093"C. In some cases, the Bodner-Partom model overpredicts the 
peak stress of the saturated loop, but the overall trend of the stress relaxa- 
tion behavior is predicted well. 

Figure 4.30 compares the model prediction and experimental results of 
stress relaxation during strain hold on the unloading portion of the saturated 
hysteresis loop of B1900+Hf at 1093°C. The Bodner-Partom calculations agree 
well with the experimental data for strain holds at . 6 % and .55%. At .5%, the 
experimental result indicates reverse stress relaxation as. the stress Is 
slightly increased during hold period while the Bodner-Partom model predicts a 
constant stress with essentially no stress relaxation. In principle, the re- 
verse stress relaxation process is capable of prediction by the equilibrium 
stress approach. However, no stress relaxation calculations are currently 
available for the Walker model. 
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FIGURE 4.28. STRESS RELAXATION DATA COMPARED WITH 
BODNER-PARTOM MODEL PREDICTION AT 
760°C. 
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Stress Relaxation at Strain Hold 



FIGURE 4.29. STRESS RELAXATION DATA COMPARED WITH BODNER-PARTOM 
MODEL PREDICTION AT 871°C. 
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FIGURE 4.30. (A) STRESS RELAXATION DATA COMPARED WITH 

BODNER-PARTOM MODEL PREDICTIONS AT 1093°C 
AND (B) STPSTS RELAXATION DATA AT THE UN- 
LOADING BRANCH OF A HYSTERESIS LOOP COM- 
PARED WITH BODNER-PARTOM MODEL PREDICTION 
AT 1093°C. 
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4.3.5 Dwell Fatigue Cycl i c Behavic r 


Comparisons of dwell fatigue data with predicted behaviors are presen- 
ted in Figures 4.31 to 4.33. A two minute hold at both the tension and com- 
pression load limits at 871 J C produces significantly more deformation than is 
predicted by either theory. Figures 4.31. Walker's model tends to produce a 
better fit of the results, although neither can be judged to be adequate, par- 
ticularly at the highest load level, +500 MPa. At 760°C, a two mi r "te hold in 
tension at 350 MPa, Figure 4.32, results in a slight amount of cyclic plasti- 
city. Walker's model overpredicts and the Bodner-Partom model underpredicts 
this quantity. At 550 MPa, both models overpredict the cyclic hysteresis 
loops. Neither the Walker nor the Bodner-Partom model is able to reproduce 
the cyclic resp-.ise with a two minute hold in compression at 982°C, as shown 
in Figure 4.33. 

4.3.6 Thermomechanical Fatigue (TMQ Cycling Constitutive Behavior 

Comparisons of predicted and experimental TMF hysteresis loops a^e 
summarized in Figures 4.34-4.36. The three the. momechan’cal fatigue (TMF) 
loading cases were run with a sinusoidally varying mechanical strain and tem- 
perature with a cyclic period of sixty seconds. The strain for three of the 
tests varied from a nominal -0.4 percent to +0.4 percent, and the temperature 
varied sinusoidally from 538°C (1000°F) to 982°C (1800°F). For the first 
loading case (Cycle A) the temperature and strain varied in phase, i.e., the 
maximum temperature occurred at the same time as the maximum strain. For the 
second loading case (Cycle B) the temperature and strain loads ware applied 
180° out of phase, and in this case the maximum temperature coincided with the 
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FtG'lRL 4.31. LOAD DWELL CYCLIC DATA COMPARED WITH THE 
WALKER AND THE BPONER-PARTOM MODEL PRE- 
DICTIONS AT 871°C. 
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FIGURE 4.32. LOAD DWELL CYCLIC DATA COMPARED WITH THE 
WALKER AND THE BODNER-PARTOM MODEL PRE- 
DICTIONS AT 760°C. 
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FIGURE 4.33. LOAD DWELL CYCLIC DATA COMPARED WITH THE WALKER 
AND THE BODNER-PARTOM MODEL PREDICTIONS AT 
982°C (A) BODNER-PARTOM, (B) WALKER. 
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FIGURE <14. THERMOMECHANICAL CYCLIC DATA COMPARED WITH THE WALKER 
AND THE BODNER-PARTOM MODEL PREDICTIONS FOR THE IN- 
PHASE TMF CYCLE. 
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FIGURE 4.35. THERMOMECHANICAL CYCLIC DATA COMPARED WITH THE WALKER 
AND THE BODNER-PARTOM MODEL PREDICTION FOR THE OUT- 
OF-PHASE IMF CYCLE. 
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FIGURE 4.36. THERMOMECHANICAL DATA COMPARED WITH THE WALKER 
AND THE BODNER-PARTOM MODEL PREDICTIONS FOR 
THE OUT-OF-PHASE TMF CYCLE WITH STRAIN HOLD: 
(A) STRESS-STRAIN RESPONSE DURING OUT-OF- 
PHASE CYCLE, (B) STRESS RELAXATION RESPONSE 
DURING STRAIN HOLD. 
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minimum strain. The third loading case (Cycle C) was a repeat of the second, 
but with the addition of a sixty second constant strain hold at -0.4 percent 
strain. 

The Walker model calculations were performed in conjunction with the 
MARC finite element deck. The structural model used was a four element axi- 
symmetric shell model using MARC element number 10, a four model quadrilateral 
ring element. The Walker model results show consistently good correlation at 
the hot temperature end of the load cycle for all three TMF missions analyzed. 
There is, however, a consistent stress overshoot at the cold (1000°F) end of 
the simulation. The Bodner-Partom model calculations, performed at SwRI, are 
in reasonably good agreement with experimental data in all aspects with the 
exception that the width of the TMF loop in Cycle C is larger than that 
observed experimentally. The stress relaxation which accompanied strain hold 
in cycle C is well predicted by the Bodner-Partom model as shown in Figure 
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5.0 HULTIAXIAL EVALUATION OF CONSTITUTIVE MODELS 


i Combined tension and torsion tests were conducted to investigate the 

effects of principal stress ratio, temperature, and loading sequence on the 
V inelastic elevated temperature multiaxial behavior of the base alloy, 

B1900+Hf. Both proportional and non-proportional loading paths were examined, 
some of which involved stress relaxation under strain hold. The various ten- 
sion/torsion experiments were designed to examine: (1) the presence or ab- 

sence of additional hardening under non-proportional (90° out-of-phase) load- 
ing paths, (2) the nature of this additional hardening, if present, (3) the 
non-axiality between the plastic strain rate and the deviatoric stress vec- 
tors, and (4) the direction of stress relaxation during strain hold. 

A considerable amount of multiaxial data have been generated in this 
program which has not yet been completely analyzed. Correlation between 
theory and experiment in this report will be directed toward examining the in- 
fluence of nonproportional loading path on the multiaxial hardening behavior 
of BI900+Hf. The results on the effects of loading "a^h on the phase angle 
between the plastic strain and deviatoric stress vectors and the direction of 
stress relaxation will be reported at a later date. 

5.1 E xperimental Procedures 

Multiaxial tension/torsion tests were performed on a servo-controlled, 
hydraulic biaxial machine with two independent feedback loops. Axial and tor- 
sional strains were independently controlled. The load capacities of the test 
facility are 44,000 N axial load and 425 N-M torque. A closeup view of the 
experimental setup is shown in Figure 5.1. Elevated temperature was provided 
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by induction heating. The biaxial test specimen is shewn in Figure 4.3 to- 
gether with the isothermal cyclic and the thermomechanical specimens. A sche- 
matic of the geometry of the biaxial test specimen can be found in Figure 4.2 

Axial and torsional strains were measured using a high temperature, 
biaxial (tension/torsion) capacitance extensometer as shown in Figure 5.2. 
The extensometer was water cooled and mounted internally to the specimen and 
grip. The coaxial roto-stator configuration of the capacitance sensing plates 
avoids the need for flexural elements in the systems and eliminates cross- 
talk. The capacitance system is of high sensitivity (6x10"^ radian and 2 ye) 
and moderate dynamic range (e < .1 sec - ^). 

The various multiaxial tension/torsion constitutive tests are 
illustrated in Figure 5.3. These tests include: (1) uniaxial loading, (2) 
proportional combined loading, (3) nonproportional loading, and (4) 
nonproportional loading with strain hold. 

Tests 1 through 4 involved incremental step strain tests in either 
torsion or tension at two temperatures, 649 and 982°F. These temperatures 
were selected to correspond to idle and takeoff temperature conditions in an 
air-cooled turbine vane. 

Tests 7 through 10 involved incremental cyclic step strain tests in 
ccnbined tension/torsion for at least two shear to extensional strain ratios 
at each temperature. The effective strain rates in these tests are approxi- 
mately 4 x 10'^ sec'l. 

Tests 5, 6, 11 and 12 were combined tension/torsion tests conducted 
under 90° out-of-phase strain cycling. The frequencies were adjusted for each 
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FIGURE 5.1. BIAXIAL SPECIMEN UNDER TEST AT 982°C. 
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(a) Rotor (upper), stator (middle) and housing (lower). 





(b) Assembled transducer. 

FIGURE 5.2. BIAXIAL TENSION/TORSION CAPACITANCE EXTENSOMETER. 
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effective strain range so that effective strain rate was approximately con- 
stant. In tests 13 and 14 the shear strain frequency doubled the axial fre- 
quency. In tests 15 and 16, the specimens were cycled 90° out-of-phase until 
saturation. A two minute strain hold was then applied at each of the strain 
peaks. After strain hold, the specimens were strained under proportional load- 
ing and the strain hold cycle was repeated. Tests 17 and 18 were similar to 
tests 15 and 16, except strain hold cycling was not applied. All the above 
tests were conducted at either 649°C or 982°C. Test 19 was a 90° out-of-phase 
nonproportional loading experiment at 871°C in which the out-of-phase loading 
was followed by in-phase loading with various strain ratio. 

5.2 Comparison of Theoretical and Experimental Results 
5.2.1 Proportional Loading Paths 

An effective stress range and an effective plastic range based on the 
van Mises criterion were used for correlating the multiaxial tension-torsion 
cyclic data. The effective stress range, Aa gff , and the effective plastic 
strain range are obtained as follow: 





where Aa and a t are the axial stress and shear stress ranges, respectively. 
Ae P and Ay P are the corresponding plastic components of the axial and shear 
strain ranges. All the in-phase multiaxial cyclic data are summarized in 
Figure 5.4 h slotting half of the effective stress range against half of the 
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effective plastic strain range. Note that all the multiaxial data obtained at 
982°C form cne small scatter band. On the other hand, the J 2 -based effective 
stress and plastic strain ranges are not able to correlate the 649°C data. 
Specifically, the shear stress range is higher than the axial stress range at 
an equivalent effective plastic stress range. The implication is that J 2 
might n )t be the appropriate criterion for B1900+Hf at 649°C. 

?oth the Bodner-Partom and Walker models are J 2 -based theories and, 
therefore, predict the same effective stress range-effective plastic strain 
range curve, regardless of the applied ay/Ae ratios. Correlations of Bodner- 
Partom model calculations and experiment, however, reveal that the overall 
agreement with the proportional biaxial data is only fair. At 982°C, the 
Bodner-Partom model overpredicts the effective stress range at a given effec- 
tive plastic strain range. The overprediction is consistent with the uniaxial 
cyclic calculations shown in Figure 4.21. The Bodner-Partom model is in fair 
agreement with the axial data at 649°C, but is rather poor with the torsional 
data. The Walker model overpredicts the plastic strain range and underpre- 
dicts the stress range so that the predicted cyclic stress strain curve falls 
below and to the right of the data at 649°C (1200°F). The agreement between 
Walker's model prediction and data at 982 g C (1800°F) is excellent. The poor 
agreement between prediction and test at 649 °C (1200 F) can be attributed in 
large part to limited test data at 649°C fr;m which model constants were 
derived (a single tensile test) and differences in material properties between 
the round bar specimens and the tube specimens used for the biaxial testing. 
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5.2.2 Won-Proportional Loading Paths 


Comparison of the multiaxial data obtained from 90° out-of-phase load- 
ing paths with those from in-phase loadings reveals no additional hardening 
under non-proportional loading paths (Figures 5.5 - 5.7). In Figure 5.5, the 
solid line represents the stress locus of the B1900+Hf specimen tested 
at = ± .4% under 90° out-of-phase until cyclic saturation. The experi- 
mental stress locus is seen to coincide with data points which corresponds to 
the stress ranges obtained from in-phase loadings at various strain paths at 
an equivalent effective strain rate and effective strain range. This finding 
of no additional hardening under nonproportional loading paths was confirmed 
by the combined nonproportional /proportional test in which the specimen was 
cycled 90° out-of-phase until saturation and then was cycled in-phase, the re- 
sulting peak stresses; o and t, lie on the same o-x locus obtained under the 
90° out-of-phase condition, as shown in Figure 5.5. The absence of nonpropor- 
tional hardening is also observed in B1900+Hf specimen cycled at = ± .6% 
at 649°C (Figure 5.6) and in specimen cycled at = ± .4% at 982°C (Figure 
5.7), and in specimen tested at 871°C. 

The Bodner-Partom model without the nonproportional hardening term, 
i.e., setting » 0 so that a = 0, was used to correlate all the 90° out-of- 
phase multiaxial data. Similarly, the Walker model without the nonpropor- 
tional hardening term Kg was used for this data correlation. As in the case 
of correlation with the in-phase data, the Bodner-Partom calculations in gen- 
eral are in good agreement with experimental data at 649°C when cycling at 
low effective strain levels (e eff = ± .4% or less), but underpredict the shear 
stress when cycling at e eff . = ± .65#. On the other hand, the Bodner-Partom 
model tends to overpredict the stresses at 982°C. The Walker model shows good 
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agreement with the axial stress at low effective strain levels (Figure 5.5), 
but underpredicts the shear stresses at 649°C (1200°F) as shown in Figures 5.5 
and 5.6. The agreement between prediction and experiment at 982°C (1800°F) Is 
quite good for Walker's model. Comparisons of the calculated and experimental 
saturated hysteresis loops of various strain ranges at 649 and 982°C are shown 
in Figures 5.8 and 5.9, respectively. The predicted hysteresis loops are gen- 
erally somewhat larger than those observed experimentally. In Figures 5.8 and 
5.9, the Walker model predictions for the saturated hystersls loops are larger 
than experiment at 649°C (1200°F), but show good agreement with experiment at 
982 # C (1800°F). 

The Bodner-Partom model predicts the frequency effect on the multi - 
axial cyclic data of B1900+Hf under 90° out-of-phase straining rather well, as 
demonstrated in Figure 5.10 and Figure 5.11. The Walker model predicts the 
frequency effect well at 982°C but somewhat overpredicts at 649°C. For non- 
proportional loading involving shear strain frequency double that for axial 
strain, the Bodner-Partom and Walker calculations also compare favorably with 
the experimental data. Figure 5.12 shows that the stresses of the calculated 
hysteresis loops for = ± .4% at 649°C are In fair agreement with the ex- 
perimental loops, but the width of the loops tend to be larger than observed. 
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FIGURE 5.10. 90° OUT-OF-PHASE MULTIAXIAL DATA AT TWO FREQUENCIES 

AT 649°C : (a) EXPERIMENT, (b) BODNER-PARTOM MODEL 

PREDICTIONS, AND (c) WALKER MODEL. 
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FIGURE 5.11. 90° OUT-OF-PHASE MULTI AXIAL DATA AT TWO FREQUENCIES 

AT 982°C: (a) EXPERIMENT, (b) BODNER-PARTOM MODEL 

AND c) WALKER MODEL. 
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FIGURE 5.12. MODEL CALCULATIONS AND NONPROPORTIONAL MULTIAXIAL DATA WITH DOUBLE FREQUENCIES. 


6.0 IMPLEMENTATION IN FINITE ELEMENT COMPUTER CODE 


The MARC nonlinear finite element computer program was the vehicle for 
Incorporating the viscoplastic models In this and two previous contract ef- 
forts (3,111. The following section will be divided Into four parts: (1) a 
brief discussion of the MARC code and the procedure for Implementing a visco- 
plastic model, (2) a review of various Integrating methods for viscoplastic 
theories In conjunction with the MARC code, (3) a summary of the conclusions, 
and (4) demonstration analysis. 

6.1 Description of the MARC Program 

In References 3 and 11, the viscoplastic constitutive theories were 
Incorporated Into the MARC program by means of an Initial stress technique. 
All of the material nonlinearity In the constitutive equations Is Incorporated 
Into an Initial load vector and treated as a pseudo body force In the finite 
element equilibrium equations. Because the viscoplastic constitutive theories 
form a "stiff" system of differential equations. It Is necessary to form the 
Incremental constitutive equation appropriate to the finite element load In- 
crement by means of a subincrement technique. 

In the subincrement technique the finite element load Increment Is 
split Into a number of equal subincrements and the viscoplastic constitutive 
theories are Integrated over the small subincrements to form an accurate re- 
presentation of the incremental constitutive equation over the finite element 
load Increment. Integration over each subincrement Is accomplished by a num- 
ber of techniques. Provided the subincrements are sufficiently small so that 
the stability level of the Integration method Is not exceeded, the technique 
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has been found to work efficiently and accurately, even for large finite ele- 
ment load Increments. However, It Is difficult for the user to pick efficient 
sublncremental steps. Thus, there Is a considerable Incentive for the devel- 
opment of an non-iterative Integration scheme which can allow large time step 
without Increasing the Instability generally associated with Intergrating the 
differential equations comprising the constitutive theory. 

The MARI code allows the user to Implement very general constitutive 
relationships Into the program by means of the user subroutine HYPELA. Within 
this subroutine the user must specify the values of the elasticity matrix D^j 
and the Inelastic stress increment vector a^ In the Incremental vector con- 
stitutive relationship. 

Ao^ s 0^j(Acj - AjoAS) - A£.j (6.1) 

The Inelastic stress Increment vector ACj Is computed In HYPELA from the vis- 
coplastic constitutive relationships. 

In Equation (6.1) a denotes the coefficient of thermal expansion and 
a j Is the vector Kronecker delta symbol. 

1 if 1 < j < 3 

( 6 . 2 ) 

0 If 3 < j < 6 

For the class of nonlinear viscoplastic constitutive relationships 
under consideration In this contract, the incremental Inelastic stress vector 
A^ depends In a highly nonlinear manner on the Incremental strain vector 
Ac ^ * 
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The solution of the incremental equilibrium is accomplished within the 
MARC code by the following algorithm. At the start of the increment the user 
subroutine HYPELA is entered to determine the elasticity matrix O^j and the 
Incremental inelastic stress vector a^.. On entry to the subroutine the input 
consists of the strain increment vector the temperature increment aT, the 
time increment At over which the incremental external load vector aP^ is ap- 
plied to the structure, and the values of the stress, strain, temperature and 
viscoplastic state variables at the beginning of the increment. Since the in- 
cremental strain vector, Ae^ = B.j au^ , can only be accurately determined 
after the solution to the incremental equilibrium relationship has yielded the 
correct incremental solution au., the strain increment vector Ae. initially 

J ' 

used U. generate A^is assumed to be the value obtained for Ae. in the preced- 
ing increment. On exit from subroutine HYPELA the elasticity matrix 0^ and 
the estimated inelastic stress increment vector Hj are obtained for each in- 
tegration point in the structure. The incremental equilibrium relationship is 
then assembled and solved for the incremental node displacement vector au^. 
The icremental strain vector, Ae. = au^., is computed and compared with 
the Initial quess for Ae. used to generate the inelastic incremental stress 
vector At j . If this incremental strain vector is equal, within a user speci- 
fied tolerance, to the incremental strain vector used to compute At. in the 
assembly phase, the solution is assumed to have converged. Otherwise, the up- 
dated strain increment vector, obtained from the solution of the equilibrium 
relations is passed into subroutine HYPELA, a new vector, Atj, is computed, 
and the equilibrium equations resolved to yield an improved value of au^. and 
Ae j. The process is repeated until the value of the vector Ae^ on the assem- 
bly phase is equal, within a user specified tolerance, to the value of the 
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vector Ae. on the solution phase. After convergence is achieved, the temper- 
ature, stress vector, strain vector and viscoplastic state variables are up- 
dated by adding the incremental values generated during the current increment 
to the values of these variables at the beginning of the increment. The pro- 
gram then passes on to the next load increment where the process is repeated. 
6.2 Integration Methods for Viscoplastic Theories in the MARC Code 

The values of D^j and at;, in the incremental constitutive relation, 

Ao. = D.j(aej - 6 j oaT) - ac. (6.3) 

are obtained by a subincrement method. Incremental values of the variables 
aT, At and (Ae. and 6 . aAT) for the current finite clement load increment are 

* J 

split NSPLIT equal values, and the constitutive equations are integrated over 
the NSPLIT subincrements to provide accurate values of D^j and ac^. Each load 
increment in a MARC analysis is divided into NSPLIT subincrements. The inte- 
gration of the constitutive equations is currently performed by using explicit 
Euler forward differences with a step size determined by dividing the MARC 
load increment by NSPLIT. 

Two methods for integrating the constitutive equations were examined: 
(1) an explicit Euler difference scheme with error estimates for revision the 
time step, and (2) integration by the NONSS (non-iterative self-adaptive solu- 
tion) method of Miller and Tanaka [12|. The forward difference integration, 
similar to [131, is based on an error estimate, E, given by 

JTaJ 2 

E = 4R + \r (6 - 4) 
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where 


sR m J§ 4t ?j se ?j 

is the inelastic strain 
S^j is the deviatoric stress 

A of a quantity is the change in the quantity over time step At. If the error 
estimate is too large 


e > ej (6.5) 

then the time step is halved and the step is repeated. If the error estimate 
is too small. 


e < £2 ( 6 . 6 ) 

then the time step is doubled for the next integration step. 

It was found to be more difficult to integrate Bodner-Partom's equa- 
tions using an explicit Euler integration and as a result the NONSS method was 
applied to the integration of Bodner-Partom's constitutive equations in the 
calculations done at SwRI. The NONSS method can be illustrated by considering 
the following two differential equations in two unknowns, 

x * f (x.y) (6.7) 
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= g (x,y) 


( 6 . 8 ) 


where ( ) = 

Consider a time step At from time t N to time t N+ j where 

At = l N+l " l N 

then the changes in each of the dependent variables are approximately 

Ax — Xj^^ — — ( il - ^) ^ x |^ + i ) (8*8) 

Ay = y N+1 - y N = At [(1 - n) y N + ny N+1 ) (6.10) 

where 0 < n < 1. 

If n is equal to zero an explicit Euler integration results while if n is 
equal to one an implicit Euler integration results. 

Equations (6.7) and (6.8) can be used for the derivatives at the new 
time steps. For example: 


Vl ■ f <*N + 1, W ’ f < X N + 4X ' y N = 4y > 


( 6 . 11 ) 



A Taylor series can be used to expand the function. Then 


X N+1 f 


* V 


+ AX f. (x 


N’ J N 


yj + Ay f,, 


( X N’ y N^ 


( 6 . 12 ) 


A similar equation holds for y N+ ^ 


» M+ 1 - 9 (*„. y M ) - 4*9. x (v y n* * W9. y <V V < 6 - 13 > 

• • • • 

[1 - nAt f, x(x N , y N )J A x - natf. Ay = At (x N - n(x N - f(x N ,- y N )]}(6.14) 

• • • • 

- nAtg, x ax + (1 - 'nAtg, (x N ,y N )]Ax = At{y N - n[y N - g (x N , y M ) 1 } (6.15) 


Equations (6.14) and (6.15) are two linear equations in the two unknown ax and 
ay. The error in the truncation of the Taylor series is included in the dif- 
ference. 


x N " ^( x N* y N^ anc * 

y N - g(x N ,y N ) 


which is applied in the next increment. 


Computer codes based on Euler integration have been completed and 
tested for both models. The computer code based on the NONSS method has not 
been completely implemented in the MARC Code. All the Bodner-Partom calcula- 
tions performed at SwRI are based on the NONSS method. 

6.3 Constitutive Model Demonstration - TMF Simulations 

The successful implementation of the constitutive models into the MARC 
computer program was established by performing uniaxial analyses and comparing 
the predicted and observed experimental response for three thermomechanical 
loading cycles. The comparison was made with both constitutive models. 

An illustrated comparison of the calculated results with the experi- 
mental data for the TMF Cycle A is shown in Figure 6.1 Comparison of Figure 
6.1 and Figure 4.35 indicates that the TMF loop based on the Bodner-Partom 
model obtained at PWA using the Euler integration scheme (Figure 6.1) is iden- 
tical to that obtained at SwRI using the NONSS method. The good agreement be- 
tween the two calculations represents a check on the computer codes developed 
at PWA and SwRI. 

6.4 Computational Efficiency 

A study to compare the stability and computing efficiency of the 
Walker and the Bodner-Partom theories was conducted. This study was done 
using the MARC finite element deck and a four element model using element type 
10, a first order quadrilateral ring element. The convergence criterion used 
internally in MARC was for the maximum residual force to be less than 0.10 
times the maximum reaction force for the model. Using this criterion, the 
finite element model was executed with three different sets of increments for 
a single TMF cycle. The computation times were recorded and are shown in 
Table 6.1. Examples of the calculated stresses are shown in Figures 6.1 and 


104 






FIGURE 6.2. COMPARISON OF EXPERIMENTAL AND CALCULATED TMF LOOPS USING 
THE MARC CODE AND 12 LOADING INCREMENTS. 






Table 6.1 


Comparison of Execution Time and Convergence Accuracy 
for Differing Increment Number per Cycle. Out of Phase TMF Cycle 


Walker 


Increments Solution Time (sec) Convergence Accuracy (X)* 


40 

11.19 

0.51 

20 

7.01 

1.49 

12 

7.41 

2.12 


Bodner-Partom 


Increments Solution Time (sec) Convergence Accuracy (%)* 



40 

35.35 

0.16 


20 

23.94 

2.43 


12 

25.04 

0.32 

* 

Estimated as 

the magnitude of the residual 

stress divided by the 


maximum stress In the cycle. 



6.2. The examples shown are tre runs using 20 Increments per cycle for both 
the Walker and Bodner-Partom theories and the runs using 1? Increments per 
cycle. The results from the 40 increments per cycle are not shown since they 
were found to be virtually identical to the cases with 20 Increments per 
cycle. The results shown in Figures 6.1 and 6.2 are a comparison of the rela- 
tive computational stability of the two theories as currently implemented in 
the MARC deck. The Bodner-Partom theory results begin to degrade at the hot 
end at 12 increments per cycle whereas the Walker theory continues to 
calculate acceptable results. Both calculations show difficulties In 
convergence whun 12 increments are used In the loading cycle, as evidenced by 
the increase in computation time per increment. This latter effect can be seen 
from the computing times shown for the 20 and 12 increnrnt cases shown in 
Table 6.1. An examination of these times shows that for both theories the 
total computing time increased slightly even though the number of loading 
increments decreased by 40 percent. 

Another significant result of this study is indicated by the run times 
shown in Table 6.1. Fo»* the same convergence criterion and the same integra- 
tion method the Bodner-Partom model required approximately 3 times as long as 
the Walker model to calculate stresses. The reason for the longer solution 
times relates to the larger number of subincrement integration steps required 
to satisfy limits on the error estimates in the HYPELA subroutine. Other 
independent studies have indicated roughly equivalent computing time for these 
and other unified models. Such comparisons will be dependent upon the 
particular problem and Integration schemes utilized. 



7.0 BENCHMARK NOTCH VERIFICATION EXPERIMENTS AND ANALYSES 


Elevated temperature testing of Instrumented notched round specimens 
was conducted to generate notch displacement data for verification of the 
analytical methodologies developed In this program. The unified constitutive 
theories of Walker and of Bodner-Partom were Incorporated Into the MARC finite 
element code for predicting structural response behavior. Calculations were 
performed based on the Walker model for the benchmark notched specimen sub- 
jected to a variety of histories and the results are discussed In this sec- 
tion. Correlation of model predictions and experimental data serves as an 
additional test of the accuracy and applicability of the constitutive models. 
The notched round specimen was selected for the following reasons: (1) to 

permit notch strains to be measured directly using conventional extensometry, 
(2) to facilitate finite element analysis by avoiding the need for 3D ele- 
ments, and (3) to verify the ability of finite element analysis to predict 
material behavior under conditions of high constraint. 

Testing was conducted for six load patterns at 871°C over load ranges 
sufficient to result In short time Inelastic behavior and over load times suf- 
ficient to Induce significant time-dependent Inelastic notch strain. The spe- 
cimen design and experimental procedures of the benchmark notch test are dis- 
cussed In Section /.l. The analytical efforts and correlation of theoretical 
calculation and experimental data are reported In Section 7.2 
7.1 Specimen Design and Experimental Procedures 

The benchmark notch specimen was designed based on the following con- 
siderations: (1) the tensile load capability of the test machine was limited 
to ± 36000 N, (2) local plastic strains needed to be measured Indirectly, to 
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within ± 20% of their cyclic value, (3) a fatigue life of 100 cycles to Initi- 
ation of a first crack was required, and (4) testing had to be performed at a 
temperature within the creep range. 

Fully reversed strain controlled testing In the NASA Host Isotropic 
Fatigue Program (NAS3-23288) has established the total and cyclic plastic 
strain levels corresponding to a fatigue life of approximately 100 cycles to 
Initiation of a first crack. At 760, 871, and 982°C the magnitude of the 
cyclic plastic strain at this life level is 0.03%, 0.15%, and 0.2%, respec- 
tively. A test temperature of 871°C was selected based on consideration that 
testing should be conduced at as high a temperature as possible to maximize 
the nonlinear effects but low enough to preclude net section plasticity. 

A notched round specimen with an elastic stress concentration factor 
of 1.5 was selected and approved for the benchmark notch experiment (Figure 
7.1). The configuration of the notched round specimen satisfied the two re- 
quirements; (1) limited net section plasticity, and (2) measurement of either 
cross notch or diametrical displacement should yield a fair measurement of 
local notch plastic strain. Cross notch and diametrical displacements were 
measured using axial and diametrical MTS extensometers. Figure 7.2 Illus- 
trates the experimental setup of the benchmark notch test. 

In order to establish the required sensitivity of the measuring de- 
vices, a nonlinear FE analyst of the specimen conf Iguratlon was performed. 
Cross notch displacement was not found to be as strong an "indication" of 
local notch plasticity as the diametrical contraction at the notch throat. A 
better measure of the notch plastic strain is the permanent radial displace- 
ment at the notch throat. For 0.15% cyclic plastic strain the extent of the 
nonlinear diametrical strain is only 70 ue. To satisfy the requirement that 
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FIGURE 7.2 BENCHMARK NOTCH ROUND SPECIMEN UNDER 








the constitutive model be able to correlate fatigue life to within a factor of 
two, the local notch strain must be measured to within i 0.03%, which trans- 
lates to a dlametrlal strain sensitivity of ± 14 ue . This level of accuracy 
is within the capability of the capacitance extensometry utilized in the 
benchmark experiment, but all other sources of error must be minimized. 

All the notched round cyclic tests were conducted at 871°C under fully 
reversed load -control led conditions at 10 cpm. The specimens were cycled at ± 
331, +352, + 365, ± 386, and ± 414 MPa for 10 cycles at each stress level. 
The load sequence was repeated twice for a total of 100 cycles. For cycling 
with dwell, the same load sequence was used with a minute hold applied at; (1) 
maximum tension, (2) maximum compression, and (3) maximum tension and compres- 
sion. In addition to the cyclic tests, a notched round specimen was pulled 
monotonlcally almost to failure at a displacement rate of 5xl0 -4 mm/sec mea- 
sured over a 25.4 mm gauge span. A monotonic tensile test with dwell was also 
conducted by step loading beginning at 379 MPa up to 517 MPa In 34.5 MPa in- 
crements. A 15 minute dwell time was applied at each load level. A summary 
of the benchmark notched round tests Is shown In Table 7.1. 

7.2 Finite Element Analyses 
7.2.1 Finite Element Mesh Selection 

The element type has been selected according to MARC recommendation. 
For the large strain plasticity problems as In monotonic loading of a 0-notch 
round specimen, simple element types (e.g., axlsymmetrlc element 10) and those 
elements containing a modified strain Integration procedure (axlsymmetrlc ele- 
ment 55) are preferable. 
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Figure 7.3 illustrates the four meshes of the test specimen considered 
in this phase of the study. The number of elements surrounding the notch 
varied from six in the coarse mesh to sixteen in the finest mesh. For pur- 
poses of mesh evaluation, a conventional time independent plasticity model was 
utilized. The assumed stress-strain material behavior reasonably approximated 
the observed tensile response of B1900+Hf at 871°C. Incremental displacements 
were applied to the model, care be^ng taken to assure that the step size was 
not excessive. Results of the four analysis were presented in Tables 7.2 
through 7.5. 

The relative computational requirements, presented in Table 7.2 indi- 
cate a ten fold increase in CPU time when the mesh density was refined by a 
factor of about two. In Table 7.3 the limit load predictions are compared. A 
difference of only 2.7 percent exists between the limit load predictions for 
the coarsest and the finest mesh. The limit load is therefore not an adequate 
criterion for determining an optimum mesh. Table 7.4 presents the predicted 
load to produce a variation in the inelastic radial displacement at the notch 
root equal to 3.05 x 10" 2 mm. As noted in the table, a difference of only 2.8 
percent exists between the load required to produce an inelastic radial dis- 
placement of 3.05 x 10“ 2 mm (1.2 x 10"^ in) for the coarsest and finest mesh. 
This parameter, a measurable quantity, is also not adequate for selecting a 
mesh. A final comparison relates the axial stress at the notch root to the 
nominal gage length strain. As noted in Figure 7.4 and Table 7.5, the mesh 
selection has a somewhat larger effect on the predicted value of the local 
stress. The intermediate mesh would appear to achieve an acceptable precision 
without incurring unacceptable CPU charges and was selected for correlation of 
the test data. 
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FIGURE 7.3. MODELS CONSIDERED IN MESH SENSITIVITY STUDY 
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TABLE 7.2 


CPU REQUIREMENTS FOR VARIOUS D-NOTCH MODELS 


| MARC Element 

Mesh 

CPU/Increment 1 
(seconds) | 

10, 4 node linear 
quadrilateral 

Coarse 

4.0 


Intermediate 

5.0 


Fine 

14.5 


Finest 

40.0 

55, 8 node quadratic 
quadrilateral 

Intermediate 

i 

23.0 


TABLE 7.3 

LIMIT LOAD COMPARISON FOR VARIOUS D-NOTCH MODELS 


| MARC Element 

Mesh 

Relative Limit 
Load 

10, 4 node linear 
quadrilateral 

Coarse 

1.027 


Intermediate 

1.0133 


Fine 

1.004 

1 1 

Finest 

1 J 

1.0 

1 _ 1 

55, 8 node quadratic 
quadrilateral 

Intermediate 

0.995 
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TABLE 7.4 






INELASTIC RADIAL DISPLACEMENT AT NOTCH ROOT COMPARISONS 
FOR VARIOUS D-NOTCH MODELS 


- , 

MARC Elemen 

Mesh 

Relative Load 
at aar of_ 
1.2 x 10‘ 3 

!# 

10, 4 node linear 

Coarse 

1.028 


quadrilateral 





Intermediate 

1.013 | 

^ 1 


| Fine 

1.003 | 

►rl- - ! 


| Finest 

1.00 | 

1 

_ _ __ _ 1 

4 

55, 8 node quadratic 

Intermediate 

0.995 

f 1 ; 

quadrilateral 




TABLE 7.5 

LOCAL NOTCH STRESS FOR VARIOUS D-NOTCH MODELS 


MARC Element 

I 

Mesh 

Local Relative 
Notch Stress 
at Limit Load 

_ 1 

1 10, 4 node linear 

| quadrilateral 

Coarse 

1.09 | 

| | Intermediate | 

1.05 | 

1 1 Fine | 

| 1.03 | 

| | Finest 

1 1 1 

1 i.o I 

1 . . 1 

55, 8 node quadratic 
quadrilateral 

Intermediate 

0.99 


i 

■ i 
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True Stress at Notch Root/6.895, MPa 
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FIGURE 7.4. NOTCH STRESS VS. GAGE LENGTH STRAIN FOR VARIOUS .VOTCH MODELS. 
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7.2.2 Experimental Results and Finite Element Predictions 

Test results, In particular, the variation of the notch root dlame- 
trlal displacement with load level, cycle number and loading wave form are de- 
scribed below. A comparison with the MARC finite element predictions using 
the Walker unified constitutive theory Is also presented. 

a) Fully Revised Cycling with No Hold (Specimen Bl) 

The specimen was cycled for ten cycles at each of five load 
levels beginning at a net section stress of i 331 MPa. No hysteresis In the 
load versus notch diametrial displacement was observed at the three lowest 
load levels. At ± 386 MPa and ± 414 MPa, some cyclic Inelastic behavior was 
noted. The magnitude of the peak to peak notch displacement and the cyclic 
notch Inelastic displace..ient did not change appreciably with Increasing cycles 
at any load level. Figure 7.5 (a) and Figure 7.5 (b). No accumulation in mean 
notch displacement with cycles was noted. Notch displacement versus load for 
the sixth cycle at + 414 MPa are presented In Figure 7.6. 

The MARC finite element code with the Walker model was used to 
simulate six cycles of loading at the ± 414 MPa load level. No hardening or 
softening with increased cycling was predicted. The peak-to-peak variation in 
notch displacement was correctly estimated; however, considerably more cyclic 
notch plasticity was predicted than observed In the test, as Illustrated in 
Figure 7.6. 

b) Fully Reversed Cycling with One Minute Hold In Tension 

(Specimen B4) 

The specimen was cycled for ten cycles at each or five load 
levels beginning at ± 331 MPa. During the first and every subsequent hold in 
tension for one minute, a notch diametrical contraction which was partially 
recovered on load reversal was noted. Notch displacement versus load response 











Specimen B1 
Throat 01 sp. Vs Load 



FIGURE 7.6. LOAD VS. NOTCH DIAMETRIAL DISPLACEMENT, TEST B1 
FULLY REVERSED. 
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4 

i 

j 

i 

for the sixth cycle at the t 331 MPa load level are presented In Figure 7.7. 
Hysteresis Is present In the results. The cyclic drift In the mean notch dla- 

; 

1 meter was small but always Increased with Increasing cycles and load level. 

The magnitude of the peak-to-peak notch displacement and the cyclic Inelastic 

1 

notch displacement. In Figures 7.8(a) and 7.8(b), shows little systematic 

i 

| variation with the cycle number. 

The finite element calculation underestimates the peak-to-peak throat 
displacement by as much as seventeen percent, but predicts the cyclic inelas- 
tic notch displacement at the ± 414 MPa load level quite accurately. No hard- 
ening or softening Is predicted. In agreement with the test data. 

The variation of the mean notch displacement at the lowest (± 331 MPa) 
and highest (i 414 MPa) load levels with hold time Is presented In Figures 
7.9(a) and 7.9(b), respectively. The notch displacement during the one minute 
tensile holds Is underpredicted by the analysis at both stress levels. 

c) Fully Reversed Cycling with One Minute Hold In Compression 




The specimen was cycled for ten cycles at each of five load 
levels with a one minute hold In compression. During the first and subsequent 
one minute compression hold, a notch dlametrial expansion which was totally 
recovered on loaa reversal was observed. This behavior is In contrast to that 
observed In Specimen 84. The notch dlametrial displacement variation with 
load for the sixth cycle at ± 331 MPa (± 48 ksl) Is shown In Figure 7.10. 
Hysteresis Is present In the load displacement response with no drift in the 
mean cyclic notch diameter at any load level. The magnitude of the peak-to- 
peak notch displacement and cyclic Inelastic notch displacement with Increas- 

i 

Ing cycles at each load level are summarized In Figures 7.11(a) and 7.11(b), 
respectively. The cyclic response stabilizes quickly upon load change. The 
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FIGURE 7.10. 


LOAD VS. NOTCH DIAMETRIAL DISPLACEMENT, TEST B5, 

FULLY REVERSED, 331 MPa, ONE MINUTE HOLD COMPRESSION. 
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FIGURE 7.11. (a) PEAK TO PEAK NOTCH DISPLACEMENT, TEST B5, FULLY 
REVERSED, ONE MINUTE HOLD COMPRESSION. 
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FIGURE 7.11. 


(b) CYCLIC NOTCH INELASTIC DISPLACEMENT, TEST B5, 
FULLY REVERSED, ONE MINUTE HOLD COMPRESSION. 
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accumulation of notch displacement with hold time and at the ± 331 MPa (t 48 
ksl) load level is presented In Figure 7.12. 

The finite element calculation overestimates the peak-to-peak 
throat displacement by as much as 44515, Figure 711(a). The cyclic inelastic 
notch displacement is also overestimated at the ± 414 MPa (± 60 ksl) load 
level by 22515, Figure 7.11(b). A constant load hold for one minute in compres- 
sion results in slightly more predicted dlametrial expansion at both load 
levels than observed in the data. A positive drift In the mean cyclic notch 
displacement with increasing test cycles and load level is predicted, contrary 
to test results in Figure 7.12. 

d) Fully Reversed Cycling with a One Hlnute hold In Tension and 

Compression (Specimen B6) 

The specimen was cycled with a one minute hold in both tension 
and compression. The variation of the peak to peak notch displacement and the 
cyclic inelastic notch displacement with the cycle number is shown in Figures 
7.13(a) and 7.13(b), respectively. No cyclic hardening is noted in the re- 
sponse. The variation with time of the notch displacement during the tensile 
hold at both ± 331 MPa (+ 48 ksi) and ± 414 MPa (± 60 ksi) load levels is pre- 
sented In Figure 7.14. Results for the compression hold are presented in 
Figure 7.15. Note that the notch strain resulting from a tensile hold was not 
totally recovered during the compression portion of the cycle and a negative 
drift in the mean hoop strain occurs at all load levels with Increasing 
cycles. 

As noted in Figures 7.13(a) and 7.13(b) the finite element anal- 
ysis predicts the peak-to-peak variation In the throat diametrial displacement 
to within 14%, but underestimates the cyclic inelastic notch displacement by 
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FIGURE 7.12. VARIATION OF NOTCH DIAMETRIAL DISPLACEMENT DURING 
COMPRESSION HOLD, TEST B5, FULLY REVERSED, ONE 
MINUTE HOLD TENSION AND COMPRESSION. 
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; FIGURE 7.14. VARIATION OF NOTCH DIAMETRIAL DISPLACEMENT DURING 

TENSILE HOLD, TEST BS, FULLY REVERSED, ONE MINUTE 
HOLD TENSION AND COMPRESSION. 
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FIGURE 7.15. VARIATION OF NOTCH DIAMETRIAL DISPLACEMENT DURING 
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nearly 60% at the highest test loads. The change in the throat displacement 
during the tensile and compressive holds. Figures 7.14 and 7.15, is well pre- 
dicted at ± 331 MPa but is underestimated at ± 414 MPa. 

e) Monotonic Tension (Specimen B3) 

The variation in the displacement at the notch throat with spe- 
cimen load level during the Initial portion of the test is presentej in Figure 
7.16. At the imposed nominal strain rate of 2xl0“ 5 sec'* the limit load was 
35.800N (8050 lb). Finite element load versus throat displacement predictions 
track the test data initially but diverge at strain conditions indicative of 
bulk yielding in the notch. The limit load is overestimated by 18%. 

f ) Stepped Toad Tension with Dwell (Specimen B7) 

The specimen was loaded to 379 MPa (55 ksi) and uploaded in 
34.5 MPa (5 ksi) increments to a maximum load of 517 MPa (75 ksi). A 15 
minute hold was applied at each load level. The change in the displacement at 
the notch throat during the initial and subsequent load changes is presented 
in Figure 7.17. The change in the displacement at the notch throat during the 
tensile hold is presented in Figure 7.18. The data is ordered well with in- 
creasing load level except at the initial load level. 

The finite element analysis performed to simulate the test 
underestimates the diametrial displacement in the notch both during the 
fifteen minute holds and during the subsequent load change by as much as 50%. 
7.3 Discussion 

Failure to correlate results of a particular benchmark notch test 
against a finite element analysis may not be attributed solely to deficiencies 
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FIGURE 7.16. LOAD VS. NOTCH DIAMETRIAL DISPLACEMENT, TEST B3, 
MONOTONIC TENSILE TEST. 
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FIGURE 7.18. PREDICTED VS. OBSERVED VARIATION OF NOTCH DIAMETRIAL 
DISPLACEMENT DURING 15 MINUTE HOLD, TEST B7, STEP 
LOAD TENSILE TEST. 
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In the form of the constitutive model, assumed model constants, or finite ele- 
ment methodology. Experimental uncertainties and bar to bar scatter In mate- 
rial properties must also be considered In the assessment. Following Is a 
discussion of these sources of error. 

The variation of the notch displacement from bar to bar at low, nearly 
elastic load levels when normalized by the apparent modulus as determined from 
the longitudinal displacement measurements, should be a reasonable measure of 
the precision with which notch strains can be measured. Uncertainty In the 
notch strains may be introduced when displacement measurements are obtained 
between two points spanning too few grains or grains of aclcular shape. The 
load-notch displacement response of each specimen was evaluated at elastic 
notch conditions and the apparent compliance established. Results, corrected 
for modulus differences, are presented in Table 7.6 for each of the six tests. 
Note that the apparent elastic diametrical compliance varies by up to 50< re- 
lative to the analytical prediction. Easy accessibility of the notch and the 
high level of reliability and precision of conventional extensometry Implies 
that the source of the variability In the notch elastic deflection may be re- 
lated to the structure of the material; in particular, the grain size, and/or 
shape of the grains in the vicinity of the notch. 

The assessment is consistent with metal lograph 1c obtained by section- 
ing of the notch specimen, shown in Figure 7.19, which reveals that the number 
of grains across the cross-section of the notch is relatively small and the 
grain size Is large compared to the notch radius. 

The scatter In creep and tensile properties from bar to bar is an Im- 
portant consideration In evaluating the finite element simulation of the test 
results. To estimate the expected variability In the benchmark test data, a 


138 



TABLE 7.6 


BENCHMARK NOTCH DATA CORRELATION WITH FINITE ELEMENT PREDICTIONS 


Sped Men 

Relative Elastic 
Compliance* 

81 Cyclic 

1.3 

B3 Monotonic 

- 

84 Cyclic 

1.49 

B5 Cyclic 

.90 

B6 Cyclic 

1.01 

B7 Monotonic 

1.37 


* Actual Value/Predicted Value 


Relative Hysteresis 
Width* 

Relative Creep 
Deflection 
During Hold* 

.20 

- 

.83 

1.27 

.62 

.76 

1.45 

2.06 

.75 

.76 


139 



ORIGINAL PAGE IS 
OF POOR QUALITY 


I 



i 

i 

i 




FICURE 7.19. A SECTIONAL VIEW OF THE NOTCH INDICATES THE PRESENCE 
OF LARGE STRAINS AT I HE NOTCH ROOT. 





statistical analysis was performed to examine the scatter in material proper- 
ties from bar to bar. Stationary strain rate data from constant strain rate 
tensile tests at different strain rates were regressed to obtain a best fit at 
871°C. A Weibull analysis of these data normalized by the regression predic- 
tions to obtain the material property scatter is shown in Figure 7.20 along 
with the the 5 and 95£ confidence limits. To verify this estimated cumulative 
distribution, an analysis was performed on ten fatigue specimens with strain 
holds at 871 C C. Of these ten, three were run with one-minute tensile strain 
holds and seven with one-minute compressive strain holds. All ten were run 
with a completely reversed strain range of .5% at a loading frequency of 10 
CPM. Two measures of inelastic material behavior were extracted from these 
tests. The measures were: 1) inelastic strain range at zero stress and 
2) the amount of stress relaxation which occurred during tne one minute con- 
stant strain hold. Each of these measures were normalized by Its average for 
the 10 tests and the results compared in Figure 7.21 with the best fit Weibull 
distribution taken from Figure 7.20. This second analysis confirms the orig- 
inal analysis of the creep and tensile data and lends assurance that the lat- 
ter may be compared to the benchmark notch data and finite element predic- 
tions. 

Table 7.6 presents two measures of the quality of the data correla- 
tion: (1) the relative width of the hysteresis loop in a cyclic test or the 
degree of nonlinear response during a load change in a monotonic test defined 
by the notch diametrial displacement versus load response, and (2) the rela- 
tive notch diametrial "creep" during the cyclic or step load test. Each entry 
represents the average value of the parameter when MARC analysis was performed 
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FIGURE 7.21. STATISTICAL ANALYSIS OF FATIGUE DATA. 
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at more than one load level. A statistical analysis of these data is presen- 
ted in Figure 7.22 and 7.23. Note that the scatter in results is consistent 
with that noted in the smooth creep and tensile data. On average, the value 
of the notch displacement during constant load holds is predicted. Figure 
7.22 ; but the cyclic inelastic displacement or total inelastic displacement 
during the monotonic tests is overestimated. Figure 7.23. Additional fully 
reversed, no hold, notched round cyclic tests should be performed to confirm 
the variability in material properties. 
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FIGURE 7.22. STATISTICAL ANALYSIS OF PREDICTED AND OBSERVED CREEP 
DISPLACEMENT AT NOTCH ROOT. 
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FIGURE 7.23. STATISTICAL ANALYSIS OF PREDICTED AND OBSERVED TOTAL 
INELASTIC DISPLACEMENT AT NOTCH ROOT. 
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8.0 HOT SECTION COMPONENT ANALYSIS 


The airfoil section of a typical cooled turbine blade has been selec- 
ted for demonstration analysis of one of the two constitutive models developed 
in this contract. Figure 8.1 depicts the finite element mesh selected for the 
analysis. The mesh includes two features which are intended to demonstrate 
the versatility of the constitutive model and its practical application to 
actual hot section components. 

Two element types are used, as shown in Figure 8.1. A segment of the 
leading edge is modelled with higher order solid elements (MARC element Type 
21) while the remainder of the airfoil is modelled with simpler elements (MARC 
element Type 7). Appropriate node tying equations are included in the finite 
element input stream to insure proper compatibility at the interface of the 
two element types. This technique of mesh refinement is commonly employed in 
structural analysis to reduce computing time and costs. Inclusion of two ele- 
ment types in the demonstration will illustrate the proper modular formulation 
of the constitutive model. 

The second feature of the model is the simulation of internal features 
of the component by modifying the density and stiffness of certain mesh elem- 
ents, Figure 8.2. The internal features being modelled in this case are small 
ribs connecting the two walls of the hollow airfoil. Of course, if detailed 
stresses and strains near these features and on a size scale comparable to the 
features are desired, a much more refined mesh would be required or the region 
would have to be subsequently analyzed using boundary conditions obtained from 
the cruder mesh. Inclusion of this feature demonstrates the ability of the 
model to be adapted for specialty use. 
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The computations using the MARC code with the Walker constitutive 
model are currently underway at PWA. The results will be included in a report 
to be issued separately. 
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FIGURE 8.2. RADIAL VIEW ILLUSTRATING INTERNAL FEATURES OF 
THE HOLLOW AIRFOIL. 


9.0 OVERALL EVALUATION OF CONSTITUTIVE MOOELS 


As Indicated earlier [1], the unified constitutive equations have some 
common properties and some essential differences. In order to evaluate their 
relative merits, it is useful to keep in mind some basic characteristics that 
such equations should possess for the purpose of engineering applications. 
Some of these are the following: 

(1) General Admissibility - the equations must be fundamen- 
tally correct from the viewpoint of the basic princi- 
ples of thermodynamics and continuum mechanics. The 
equations should be consistent with the fundamental 
physics of elastic and inelastic deformation mechanisms 
and thereby be able to represent response properties 
over a wide range of strain rates, temperatures, and 
loading conditions with minimum alteration of the' basic 
equations. 

(2) Range of Applicability - the equations should be cap- 
able of representing essential material response char- 
acteristics such as strain rate sensitivity of plastic 
Flow, microplasticity ( 1 . e . , small plastic straining at 
low stress levels), isotropic and directional harden- 
ing, thermal recovery of hardening, creep and stress 
relaxation, cyclic hardening and softening, cyclic 
creep and relaxation, temperature and pressure depen- 
dence of plastic flow, and, if applicable, temperature 
history effects. 
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(3) Material Data Requirement - the equations should con- 
tain a relatively small number of Identifiable material 
constants which are determinable from standard test 
data. It Is preferable, but not essential, that the 
material constants be amenable to physical interpreta- 
tion. 

(4) Accuracy - the equations must be capable of making ac- 
curate predictions of constitutive behavior. A high 
level of accuracy is required not only for the simple 
cases to which the constitutive models are fitted out 
also for cases of complex loading histories for which 
data are not used in the determination of the material 
constants. 

(5) Integrability - the equations should be capable of 
being incorporated into standard finite element struc- 
tural programs in a direct manner. 

These five cr* eria will be used to evaluate the overall performance 
of the Bodner-Partom and the Walker models. 

(1) General Admissibility : 

Both the Bodner-Partom and the Walker mode's appear to be funda- 
mentally correct from the mechanics and thermodynamics viewpoint and are con- 
sistent with the overall elastic and inelastic deformation mechanisms. The 
modified Walker model utilizes an exponential function for the kinetic equa- 
tion relating the inelastic strain rate to the stress and the internal vari- 
ables. This function is similar, but not Identical, to that used in the 
Bodner-Partom model. The form of the evolution equation for the isotropic 
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hardening variable Is also similar for both models, with the exception that 
thermal recovery Is not included In the Walker model. The main differences 1*i 
the two models are In the flow law and the treatment of directional hardening. 
Comparison of model calculations with experimental data, however, reveals no 
practical difference In the ability of both models to predict cyclic behavior 
of B1900+Hf subjected to uniaxial and multlaxlal loading under both propor- 
tional and nonproportional strain paths. 

Theoretically, the Walker model, which is based on the equili- 
brium stress, is expected to be capable of modeling reversed creep, reversed 
stress relaxation, and a phase angle between the Inelastic strain rate and the 
deviatorlc stress vectors. On the other hand, stress relaxation and nonpro- 
portional multiaxlal calculations and test results indicate that the lack of 
these potential capabilities in the Bodner-Partom model are not severe liabil- 
ities. For example, the stress relaxation results of the Bodner-Partom model 
are In quite good agreement with the experimental data for B1900+Hf even 
though the small amount of reversed stress relaxation observed In the data is 
not predicted. The practical importance of the phase lag in non-proportional 
loading requires further evaluation. 

(2) Range of Applicability : 

On the basis of this investigation, it appears that the range of 
applicability of the Bodner-Partom and the Walker models Is quite broad. It 
is noted that the present Investigation is limited to the relatively slow to 
moderate strain rate range (10' 7 - 10'^ sec"*) and that extension of the study 
to high strain rates would be of interest. For uniaxial loading, both models 
can be used for modeling: (1) monotonic stress-strain behavior, (2) creep be- 
havior, (3) cyclic stress-strain behavior, (4) rtress-relaxation behavior, and 
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(5) thermomechanical constitutive behavior. For multlaxlal loading, both 
theories are capable of modeling Inelastic response under complex loading his- 
tories involving either proportional and/or nonproportional loading paths. 
Because of this wide range of applicability, both theories demonstrate that 
the unified concept Is a viable approach for moiling time-temperature depen- 
dent inelastic deformation. The Walker model, howover, suffers from the draw- 
back that it requires considerable modifications ..ru.n the theory Is applied to 
different structural alloys. This deficiency is evidenced by comparing the 
version of the Walker model for Hastelloy-X [3| to that for the B1900+Hf usrd 
in the present program. 

(3) Material Data Requirement : 

Both the Bodner-Partom and the Walker theories require only a 
limited data base for * jdel constant determination. All the Bodner-Partom 
model constants, except a^, can be evaluated from uniaxial tensile data. The 
multlaxlal data of B1900+Hf Indicate no additional cyclic hardening due to 
nonproportional loading. Since alditional cyclic hardening resulting from 
nonproportional loading has been observed in Hastelloy-X ! and other struc- 
tural alloys [14,15|, the ncnproportional hardening effect i:> therefore mate- 
rial dependent and must be evaluated on an individual basis. Thus, it appears 
that some nonproportional multiaxial data would be required for evaluation of 
certain material constants in any unified constitutive model, e.g. in the 
Bodner-Partom model which equals to zero in this case. Material constants in 
the Walker model can be evaluated from uniaxial tensile, creep and cyclic 
data, while some nonprcportional multlaxlal data would also be required for 
including the nonproportional hardening effect. 



The number of material constants that needs to be evaluated from 
experimental data is roughly the same in the Bodrier-Partom and Walker models. 
A difference between these two theories is that most of the materials con- 
stants in the Bodner-Partom theory are temperature- independent, while most of 
those in the Walker theory are temperature-dependent for B1900+Hf. In addi- 
tion, the Bodner-Partom mode! constants can be evaluated more readily and at- 
tributed to specific mechanisms. The consequences are: (1) the procedures 
for evaluating the Bodner-Partom model constants are simpler than for the 
Walker model, and (2) the data requirement is also less for the Bodner-Partom 
model. 

(4) Model Accuracy : 

Despite the fact that only tensile data were used for evaluating 
material constants, the Bodner-Partom model is as accurate, if not more accu- 
rate, than the Walker theory. Both models performed equally well in the ten- 
sile data correlation, with the difference that strain aging is included in 
the Walker model but not in the Bodner-Partom model. In principle, a strain 
aging term c n also be included in the latter model. This high level of ac- 
curacy is expected because most of the material constants were determined 
using the tensile data. The Bodner-Partom model appears to be more accurate 
than the Walker model in correlating the steady state creep behavior of 
B1900+Hf. Fur correlating cyclic response, both models predict the stress and 
the plastic strain ranges equally well. The shape of the hysteresis loops and 
the drift of the mean stress are better predicted by the Bodner-Partom model. 
A rigorous test of the two models appear to lie in .he TMF cycle. In all the 
three TMF cycles examined, both models predict the shape, the stress range, 
and the stress relaxat.on behavior associated with the TMF hysteresis loops 
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quite well. Note that the Bodner-Partom model Is quite effective in predict- 
ing the stress relaxation behavior of B1900+Hf under both isothermal and TMF 
conditions. The success of these models in predicting the TMF hysteresis 
loops suggests that no additional terms need to be included in the unified 
constitutive models specifically for the TMF histories. Only the static ther- 
mal recovery terms and the variation of material constants with temperature 
were used in predicting the TMF response of B1900+Hf. Further study of non- 
isothermal behavior is included in future work. 

Both models are able tn predict the first cycle mean stress of 
the stable hysteresis loops of B1900+Hf cycled under an imposed mean strain. 
In the case of the Bodner-Partom model, the mean stress at the sixth cycle is 
also predicted in most instances, but the model appears to have difficulties 
of maintaining a mean stress at small inelastic strain ranges. Under these 
conditions, both models have difficulties in predicting the drift of the small 
mean stress over a large number of cycles (10-20 cycles). 

Neither model is capable of predicting the hysteresis loops with 
dwell under load-controlled conditions. The discrepancies can be explained on 
the basis that load-controlled tests are more sensitive to material variation. 
Hardening of B1900+Hf generally saturates rather rapidly and the stress-strain 
curves level off after a few tenths of one percent of strain so that a slight 
variation in the material properties would cause a large difference in the 
plastic strain under load-controlled conditions. As a result, the unified 
models have difficulties in predicting creep strain as function of time in 
both the primary and steady-state creep regimes, even though the steady-state 
creep rate, which is the slope of the creep strain vs. time curve, is pre- 
dicted fa’rly accurately. 



The viscoplastic model analysis did not yield an adequate pre- 
diction of the benchmark notch behavior of B1900+Hf. Material variability is 
certainly one aspect of discrepancies which need to be considered. Model ac- 
curacy may also be a factor. It appears there is a need to perform the same 
benchmark notch calculations using the Bodner-Partom model in order to have a 
one-to-one comparison and to examine whether better correlations could be ob- 
tained. 

(5) Integrability : 

The TMF calculations indicate that the Bodner-Partom model re- 
quired approximately 3 times as much as computer time as the Walker model for 
performing the same calculations. The stiffness of the Bodner-Partom model is 
generally attributed to the exponential kinetic equation. Since the modified 
Walker model also uses a similar, though not identical, exponential equation, 
it is not clear why the Bodner-Partom calculations require more CPU time. In 
addition, the Walker model contains more internal variables and evolution 
equations. The issue of the integrability of the unified models needs to be 
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10.0 SUWARY AND CONCLUSIONS 


10.1 Summary 

During the base program, the following tasks were accomplished and re- 
sults obtained: 

o A literature review of existing unified constitutive 
models for time and temperature dependent inelastic be- 
havior of high temperature alloys was completed. All 
models contain three essential components - a flow law; 
a kinetic relationship between stress, inelastic strain 
rate, internal variables, and temperature; and evolu- 
tion equations for the hardening variables. Ten dif- 
ferent models were compared relative to these compon- 
ents. 

o A review of numerical integration procedures applicable 
to these "stiff" equations was completed. 

o The particular models of Bodner and Partom and of 

Walker were selected for detailed study and implementa- 
tion in a finite element code. 

o The alloy selected for testing under th ase program 

was cast B1900+Hf. Specimens were fabricated at PWA 
from a single heat of this material, which included 
tensile, creep, cyclic, TMF, biaxial, and benchmark 
notch round specimens. 
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Mechanical tests including tensile, creep, cyclic, 
stress relaxation, TMF and biaxial tests were conducted 
to generate uniaxial and multiaxial data for correla- 
tion and verification of the unified constitutive 
models developed in this program. 

A systematic approach was developed for evaluating 
material constants for the Bodner-Partom and for the 
Walker models. 

Both the Bodner-Partom and the Walker models were used 
to correlate the uniaxial and multiaxial behavior of 
B1900+Hf. Correlations include tensile, creep, cyclic, 
stress relaxation, TMF, and biaxial deformation under 
simple and/or complex histories. Predictions of the 
two models are generally in good agreement with experi- 
mental data. 

Benchmark notch round data were generated and corre- 
lated with the Walker model. Agreement between theory 
and experiment is less than desirable. Material pro- 
perty variation may be a factor in influencing the 
agreement between theory and experiment. 

Both models were implemented in the MARC finite element 
computer code. The in-phase TMF cycle was calculated 
using both the Bodner-Partom and the Walker models by 
means of an Euler integration scheme. The former model 
was found to be numerically less efficient than the 
latter model in this exercise. 
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A hot section component is being analyzed using the 
Walker model in conjunction with the MARC finite- 
element computer code. The component being analyzed is 
the airfoil section of a typical cooled turbine blade. 

A separate report will be issued describing this anal- 
ysis. 

10.2 Conclusions 

The results of this program provide strong evidence for the applica- 
bility of the unified constitutive equation approach to describe the strongly 
nonlinear, time and temperature dependent response of metals to arbitrary load 
or deformation histories. We believe that the current effort, perhaps the 
most comprehensive verification of plasticity models to date, has demonstrated 
positively the essential feature of the "unified" model, i.e., that all as- 
pects of inelastic deformation can be incorporated in a single term. This 
term is based upon a tensorial relation (flow law) between the deviatoric com- 
ponents of inelastic strain-rate and stress, a scalar kinetic relation ex- 
pressing the dependence of deformation rate on stress, internal variables and 
temperature and a set of evolution equations describing both isotropic and 
directional strain or work hardening. As a minimum, the two models studied 
were correlative in that they demonstrated reasonable agreement with all the 
experimental data generated over a very wide range of temperature, deformation 
rate and load cycle. This was accomplished with a fixed model and fixed set 
of material constants. In a broader sense, the models were predictive in that 


the functional representations and the material constants were obtained from a 



small subset of the total tests performed. Even better overall correlation 
could have been achieved If all the data were used to optimize the functions 
and constant: in each model. Some more detailed conclusions follow. 

A major criticism of the unified models has been the apparent (and 
often real) difficulty in defining and experimentally determining the multiple 
material constants employed in the models. We believe we have made signifi- 
cant progress in this regard. The different terms in the constitutive equa- 
tions are very interactive so that it is essential to have a physical insight 
into the meaning of each function and its associated constants. Procedures 
presented in this report delineate a deterministic method for determining the 
constitutive constants in a sequential manner from a relatively small material 
test data base. In the Bodner-Partom model, which is functionally less com- 
plex and has fewer material constants than the Walker model, it was possible 
to Identify all the constants for B1900+Hf from a set of monotonic tensile 
stress-strain curves. Future work will determine if this procedure is more 
generally applicable and if it is possible to automate the constant determina- 
tion procedure for a given model. 

The additional complexity of the Walker model (and other functionally 
similar models) resides primarily in the inclusion of the "back stress" or 
more appropriately designated "equilibrium stress" tensor in the flow law. 
This is a controversial area. The Bodner-Partom formulation, which does not 
employ the equilibrium stress formulation, has been shown in this study to be 
equally satisfactory (in some cases better) than the equilibrium stress form 
employed in the Walker model. This parity is achieved with fewer material 
constants. However, more analysis of the data and study is needed of the vec- 
torial relationship between the incremental strains (or strain rate) and the 



deviatoric stress. These relationships will be more sensitive to the detailed 
flow law than are the magnitudes of the independent variables. The biaxial 
data for B1900+Hf also show some indication of departure from the J 2 based 
flow criterion assumed in both models. A more complex j£, J 3 criterion would 
additionally affect the directional aspects of incremental strains with re- 
spect to the applied stress state. 

The question as to whether the time integral of plastic strain rate or 
of plastic work rate or of some other deformation related parameter is the ap- 
propriate measure of current state (hardening) is left open. It may well be a 
material dependent parameter. However, the choice does result in some spe- 
cific differences in predicted response such as in cross hardening and de- 
formation history dependence. 

The current results indicate that properties determined from isother- 
mal tests yield quite accurate results for non-isothermal strain cycling 
tests, even when the temperature cyclic range includes thermal recovery pro- 
cesses. This result, if general, would reduce the material testing require- 
ments to isothermal only. The need for biaxial testing at this point appears 
to be material dependent, some alloys showing significant additional hardening 
under non-proportional biaxial load cycling while B1900+Hf demonstrated equi- 
valent hardening under proportional and non-proportional cycling under the 
conditions tested. Additional work is needed to determine whether a general- 
ized model can be developed based upon uniaxial test data only. Biaxial test- 
ing significantly increases the data generation complexity and cost and, 
therefore, is not routinely available. Fortunately, the non-proportional bi- 
axial effects are small In most practical applications so that the uniaxial 
data provides a close approximation. 


While the correlative and predictive capacity of the unified models 
has been demonstrated, their practical application will depend on their adapt- 
ability in efficient computational algorithms when incorporated into finite 
element or other numerical methods for structural analysis. Preliminary use 
with the MARC finite element code indicates that some improvement is desirable 
in the numerical integration procedures used with these stiff, nonlinear equa- 
tions. The computation time for an equivalent problem is approximately the 
same order of magnitude as when using classical, time-independent plasticity 
formulations. 
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APPENDIX 


The main difference in the various unified theories is th r treatment 
of directional or kinematic hardening (ij. Differences exist not only in the 
choice of the flow law but also 1r the evolution equations of the internal 
variables. Two flow laws which are commonly used in the unified models are: 


and e i j = x 2 < S ij " n ij ) 

with = 0; S^. = - ~y- «ij ° kk 


where x^ and x^ are deformation-history-dependent parameters. and . are 
the deviatoric and direct stresses, respectively. The tensor represents 
the "equilibrium stress," the directional hardening internal variable which 
has also been referred to as the "back stress" and the "rest stress." 

For uniaxial stress, the flow law represented by equation (I) can be 
rearranged as: 


o = f ] (c P , T) (K ♦ D) 


( 3 ) 


where K and 0 are the isotropic and directional nardening internal variables, 
respectively for this form of the flow law. The history dependent parameter, 


f^, Is a function of Inelastic strain rate and ♦■emperature. In a similar 
manner, the equilibrium-stress-based flow law, Equation (Z), can be expressed 
as: 


o * f 2 (e P , T) K + Q (4) 

or -°-j^ = f 2 (e P , T) (5) 

where f 2 Is also a function of Inelastic strain rate and temperature and Is 
generally different from f^. An interesting observation of Equation (3) Is 
that stress and Internal variables can be scaled using the function .'j. On 
the other hand. Equation (4) Indicates that stress cannot generally be scaled 
using the function f 2 ; instead, it is the ratio of the ((a~a)/K| which Is 
scaled by f 2 * Thus, the manner in which the direction hardening variable is 
Incorporated in the flow law appears to have significant ramif ication* on con- 
stitutive behavior. The effects are expected to be more complex at elevated 
temperature because of concurrent hardening and thermal recovery. 

The objective of this appendix is to present a comparison of the dif- 
ferent approaches for modeling directional hardening Including thermal recov- 
ery effects. In addition to the flow law, the choice of the hardening measure 
will also be examined; the plastic work rate will be used In conjunction with 
Equation (3), whereas the plastic strain rate will be used with the equili- 
brium stress approach. Equation (4). A procedure for delineating isotropic 
and directional hardening from uniaxial data will be given. 
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(A) Low and Interned late Temperature Behavior without Thermal Recovery 

The evolution equations for tne Isotropic and t*e directional harden- 
ing Internal variables without thermal recovery are generally written In a 
saturation-type format, which can be viewed as a c^hination of a linear hard- 
ening and a "dynamic recovery" term: 


* • 

K = m 1 (Kj - *:)M (6) 

• « 

q * m, (Qj - n)M , or, alternatively, (7) 

• • 

0 = m 2 (0 1 - D)M 

where and are material constants. Kj and flj (or Oj) are the limiting 
values for the Isotropic and the directional internal variables, respectively. 
The hardening measure, M, can be either the 1**V.as;nc strain .ate or the in- 
elastic work rate. The dot above the Internal variables denotes the rate of 

charge with time. 

The approach taken in Equation (3) using the plastic work rate as the 
hardening measure Is first considered. A particular unified model which util- 
izes this approach is the Bodner-Partom theory (2|. In the absence of thermal 
recovery and on the physical bas<s that 0 saturates much more rapidly than K 
for most materials, l.e., m ? > m^, the following approximations to Equation 
(3) can be made In the region of small plastic strains. 
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K = (o/f x ) - 0 . (o/fj) - 0j 

(8) 

and, 

D = (c/f x ) - K * (a/f.) - K q 

OO 


where K Q is the initial value of K. Substituting these approximations in the 
general expression for y = da/dW p obtained from the evolution equations and 
Equation (3), namely. 


Y = f^flCj-K) + m 2 (0 t - D) ] (10) 

leads to 

y = f| !m^(l\j + 0^ + m^ (K q + D^) j - (m^ + m 2 ) a (11) 

This indicates that a plot of y against a should be linear with a slope of 
+ m 2 in the range of small plastic strains. 

A t-he larger strains, D fully saturates to the limiting value so 
th t setting D = and again using Equation (d) for K gives 

y = f 1 m 1 (K 1 * D^) - mjo (12) 

which is a linear relation with slop . These relations indicate that the 
y-o function should be essentially .inear at both the low and high plastic 
strain regimes. 

For models based on the equilibrium stress approach and the plastic 
strain rate as the hardening measure, the general expression for e = dc/de p 
obtained from the evclutior equations and Equation (4) is: 
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f 2 - K ) + - n) 


(13) 
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Using Equation (A) t K and a at small plastic strains can be approximated as: 


K = (o - a)/f 2 * (o - a 1 )/f 2 


(14) 


a = a - f,,K = a - f ~K 
c Z 0 


(15) 


Substituting these approximations into Equation (13) leads to 


e ■ (f 2 K 1 + Qj) + m 2 (f 2 K q + a^ - (n^ + m ? ) a (16) 


Assuming a saturates to a^ at large plastic strains, setting a = flj and sub- 
stituting Equation (14) for K into Equation (13) leads to 

0 = m 1 (f 2 Kj + a^ - rr^o (17) 

Equations (16) and (17) indicate that the 0-0 C'rve should be linear with a 
slope of tn^ + m 2 at low plastic strains ?nd with slope m^ at high plastic 
strain. 
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I (B) High Temperature Behavior with Thermal Recovery 

| The effect of static thermal recovery can be included by incorporating 

thermal recovery terms in the evaluation equations of the Internal variables. 
The Orowan-Bailey hardening/recovery format is generally used (11; 
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K = m x (Kj - K) M - Aj Kj 


K - K, 


V 1 


,0 , r 2 


(18) 


D = m 2 (D 1 - D) M - A 2 Kj |J-| ‘ (19) 

where Kg is the fully recovered value of K at a particular temperature T, and 
l<2 can generally be taken to be K Q , the initial value of K. 

The general expression for y including thermal recovery can be ob- 
tained by differentiating Equation (3) with respect to Wp and substituting 
• • 

for K and 0 , from Equations (18) and (19), giving 
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For small plastic strains, the approximation Eauations (8) and (9) are substi- 
tuted into tquat on (20) and app< jximating Dj * 0 = (o-o 0 )/fj in the recovery 
term, leads to 
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with cj = f , K (T) 
0 1 0 


( 22 ) 


Assuming A = A^ = Ag and r = r^ = rg the rate of work hardening including 
thermal recovery becomes: 
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Usirg a similar procedure, the rate of strain hardening for the equilibrium- 
stress-based models using the plastic strain rate measure can be derived: 


0 = m 1 (f 2 K 1 + flj) + m^ (f 2 K q + Qj) - + m 2 ) o - 


A (f 2 K 1 ) 1 ' r 

|P 



(24) 


It is noted that when A is zero. Equations (23) and (24) reduce to Equations 
(11) and (16), respectively. This occurs when static thermal recovery is un- 
important. Furthermore, Equations (23) and (24) indicate that the slopes of 
the linear portions of either the y-o or the e-o curves are drastically 
changed when thermal recovery is prominent. Thus, it is possible to identify 
the temperature regime at which static thermal recovery is important by 
comparing the y-a or e-o data over a temperature range. 

An examination nf Equations (11), (16), (23), and (24) also reveals 
that the hardening expressions derived for the modeling approaches with or 
without the equilibrium stress concept are of similar forms. In both cases, 
the hardening rates, y and 9, decrease with increasing stress linearly at the 
small and large plastic regimes when thermal recovery is not important and the 
slopes of these curves are influenced by thermal recovery. Significant 
results from the theoretical analysis a.e the observations: (1) the scaling 
of a and (K+0) by the function f^ in the approach without the equilibrium 
stress and the scaling of (o-n)/K by another function f 2 in the equilibrium 


stress approach; ( 2 ) scaling of the y - 0 or 9-0 curves by strain rate and tem- 
perature through their influence on f^ and f 2 , and (3) the possibility of de- 
lineating the isotropic and directional hardening components from uniaxial 
hardening data on the basis of the linear portions of the y-a or e-o. It 
should, however, be noted that despite many similarities, the two hardening 
approaches are different and their respective appli. ability should be examined 
with aspect to the experimental data. Comparison of the y-o and 0-0 curves 
indicates that the former formulation is more appropriate for B1900+Hf, as 
demonstrated in Figure 4.7. 


